Primordial Perturbations in Einstein- Aether and BPSH Theories 
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We study the primordial perturbations generated during a stage of single-field inflation in Einstein- 
aether theories. Quantum fluctuations of the inflaton and aether fields seed long wavelength adi- 
abatic and isocurvature scalar perturbations, as well as transverse vector perturbations. Geomet- 
rically, the isocurvature mode is the potential for the velocity field of the aether with respect to 
matter. For a certain range of parameters, this mode may lead to a sizable random velocity of the 
aether within the observable universe. The adiabatic mode corresponds to curvature perturbations 
of co-moving slices (where matter is at rest). In contrast with the standard case, it has a non- 
vanishing anisotropic stress on large scales. Scalar and vector perturbations may leave significant 
imprints on the cosmic microwave background. We calculate their primordial spectra, analyze their 
contributions to the temperature anisotropies, and formulate some of the phenomenological con- 
straints that follow from observations. These may be used to further tighten the existing limits on 
the parameters for this class of theories. The results for the scalar sector also apply to the extension 
of Hofava gravity recently proposed by Bias, Pujolas and Sibiryakov. 



I. INTRODUCTION 

The enigmas of dark matter and cosmic acceleration 
have motivated the exploration of theories where grav- 
ity is "modified" at large distances. On the other hand, 
the range of possibilities for constructing such theories is 
severely limited by the requirement of general covariance, 
and for that reason, most of the proposed alternatives to 
General Relativity (GR) can in fact be cast as GR cou- 
pled to new fields. 1 

Cosmic acceleration may be due to a scalar field slowly 
rolling down a potential [H [5] , or simply sitting in one of 
its local minima [5] . Alternatively, it can be driven by the 
non-minimal kinetic term of a k-essence scalar field with 
a Lagrangian of the form p(X,(f>), where X = d^d^cj). 
This form is quite versatile, and can be used to mimic 
cosmic fluids with a wide range of possibilities for the 
effective equation of state and speed of sound, including 
those which are characteristic of dark energy and cold 
dark matter [ID] . 

The gradient of the k-essence field, is a time-like 
vector which spontaneously breaks Lorentz invariance, in 
a way that is parametrically independent of its effects on 
the time evolution of the background geometry. In par- 
ticular, Lorentz invariance can be spontaneously broken 
by dfj,(j) while the background spacetime remains maxi- 
mally symmetric, a situation which is known as ghost 
condensation [IT] . Still, the "fluid" responds to the grav- 
itational pull of ordinary matter, leading to modifications 



A counterexample is the DGP brane-world scenario, where grav- 
ity is modified in the infrared by a continuum of Kaluza-Klcin 
gravitons pp. Because of the continuum, DGP cannot be formu- 
lated as a standard four dimensional GR with additional fields. 
See also [2] and [3], for recent related proposals in the four di- 
mensional context. 



of the long range potentials. 

More generally, theories with a massive graviton can 
be written in a covariant form as GR coupled to a set 
of "Stiickelberg" scalar fields 4> A with non-minimal ki- 
netic terms, whose gradients have non-vanishing expec- 
tation values [TH [12] • Depending on the interactions and 
the expectation values of the condensates, this can de- 
scribe different phases of massive gravity. Aside from the 
Lorentz preserving Fierz-Pauli phase [12] (see also [14] ). 
Lorentz breaking phases have been investigated in [13] . 
Some of these have interesting phenomenology, such as 
the absence of ghosts in the linearized spectrum, a mas- 
sive graviton with just two transverse polarizations, and 
weak gravitational potentials which differ from those in 
standard GR by terms proportional to the square of the 
graviton mass [TJ2 EH US] • 

Additional fields of spin 2 have been considered in 
bi-gravity (or multi-gravity) theories 117] , where space- 
time is endowed with several metrics interacting with 
each other non-derivatively. Due to general covari- 
ance, only one of the gravitons in the linearized spec- 
trum stays massless, while the remaining ones acquire 
masses proportional to the non-derivative interaction 
terms. Lorentz invariance can be broken spontaneously 
even in cases where all metrics are flat, provided that 
their light-cones have different limiting speeds. This 
leads to phenomenology [IS] similar to that of certain 
phases of Lorentz breaking massive gravity referred to 
above [131 E3 HE] ; of which multigravity can be seen as 
a particular realization. 

Finally, additional vector fields have received consid- 
erable attention in cosmology. Effective field theories 
for vectors are strongly constrained by stability require- 
ments. Typically, those with non-trivial cosmological dy- 
namics contain a massive ghost 19J, which can be re- 
moved from the spectrum by sending its mass to infinity. 
This amounts to imposing a fixed-norm constraint on the 
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vector, which in turn forces a Lorentz-breaking vacuum 
expectation value. This led Jacobson and Mattingly to 
dub this type of models Einstein-aether theories [2U1 HI] • 
Their low-energy excitations are the Goldstone bosons of 
the broken Lorentz symmetry, 2 which will participate in 
the dynamics of the long range gravitational interactions. 

An interesting recent development is the proposal by 
Hofava |25j that a Lorentz-breaking theory of gravity 
may be renormalizable and UV complete. The break- 
ing of Lorentz invariance in this case is implemented by 
introducing a preferred foliation of space-time, but no 
additional structure. As pointed out in [26], any the- 
ory with a preferred foliation can be written in a gen- 
erally covariant form by treating the time parameter 
which labels the different hypersurfaces as a Stiickelberg 
scalar field T. The foliation is considered to be physical, 
but not the parametrization, and therefore the covari- 
ant theory should be invariant under field redefinitions 
T — > /(T). In other words, the Lagrangian can have a 
dependence on the unit normal to the hypersurfaces, but 
not on the magnitude of the gradient T,n (in contrast 
with the examples of k-essence and ghost condensation 
mentioned above). From this observation, Bias, Pujolas 
and Sibiryakov showed 27J that Hofava gravity could be 
extended by including in the action all terms compatible 
with reparametrization symmetry, and consistent with 
power counting renormalizability. Interestingly enough, 
this extension also cured certain problems in the scalar 
sector of the original proposal (such as instabilities and 
strong coupling at low energies [33]). Jacobson [25], has 
recently clarified the relation between the Einstein-aether 
theory and this extended version of Hofava gravity, which 
he dubbed BPSH gravity. In particular, he pointed out 
that any solution of Einstein-aether where the vector field 
is hypersurface orthogonal is also a solution of the low 
energy limit of BPSH gravity. 

Since the aether only interacts gravitationally, any sig- 
nal of it must be proportional to a power of (E/Mp) , 
where Mp is the reduced Planck mass, and E is an en- 
ergy scale. Thus, even though the aether contains mass- 
less fields, its presence is hard to detect. In that re- 
spect, inflation provides an interesting window to probe 
the aether and its implications. During inflation, short- 
scale vacuum fluctuations of light fields are transferred 
to cosmological distances, where they may leave an ob- 
servable imprint. It is thus natural to look for signatures 
of Einstein-aether on the spectrum of primordial pertur- 
bations, which is the subject to which we devote this 
article. 

Previous work on this subject has been done in Refs. 
[29), 130] , although in a somewhat narrower region of pa- 



2 In theories with spontaneously broken spacetime symmetries, the 
number of Goldstone bosons does not generally agree with the 
number of broken generators. However, if the order parameter 
that breaks the spacetime symmetry is spacetime-independent 
(as the constant aether field), then both numbers do agree |24| . 



rameter space and with somewhat different conclusions. 
In the scalar sector, we find that there is a primordial 
isocurvature mode, which can be interpreted as the ve- 
locity potential for the aether with respect to matter. De- 
pending on the aether parameters, this mode can grow 
on superhorizon scales, leading to a large random ve- 
locity field for the aether. Similar results apply to the 
transverse vector sector. These perturbations may thus 
be of phenomenological interest. We also find that the 
isocurvature mode is strongly correlated with the usual 
adiabatic mode, which corresponds to curvature pertur- 
bations in the co-moving slicing. 

For previous work on the impact of adiabatic scalar 
perturbations on the cosmic microwave background ra- 
diation (CMB) and large scale structure in (generalized) 
aether theories, see [5T1 13"2] . 

While this paper was being prepared, an interesting re- 
lated paper by Kobayashi, Urakawa and Yamaguchi ap- 
peared [33] . which analyzes the post-inflationary evolu- 
tion of the adiabatic scalar mode in BPSH theory. Where 
we overlap, our conclusions agree. 

The plan of the paper is the following. In Section [TT| we 
review the basics of Einstein-aether theory and the ho- 
mogeneous cosmological solutions. Sections |III[ |TV| and 
|V| are devoted to the analysis of tensor, scalar and trans- 
verse vector perturbations respectively. Section |VI| ana- 
lyzes the contribution of vector modes to the CMB spec- 
trum. 

Readers familiar with the details Einstein-aether or 
BPSH theory are encouraged to jump directly to the con- 
cluding Section [VII| for a self contained summary of the 
main results. 

Appendix A summarizes the existing bounds on the 
parameters of Einstein-aether theories. Appendix B dis- 
cusses the equations of motion for the scalar sector of 
the theory in the longitudinal gauge. Appendix C deals 
with the canonical reduction of the scalar sector to the 
two physical degrees of freedom (a necessary step for the 
proper normalization of the vacuum fluctuations). Ap- 
pendix D contains a derivation of the long wavelength 
adiabatic and isocurvature scalar modes, for generic mat- 
ter content and expansion history. Appendix E derives 
the CMB temperature anisotropies due to vector modes. 

II. EINSTEIN-AETHER THEORIES 

The Einstein-aether is described by the most general 
Lagrangian with two derivatives acting on a vector field 
of constrained norm [20] . 

C A = c 1 V a A' f V a A~ / +c 2 V a A a V~ / A'»+c 3 V a A~ t V~ / A a - 
-c i A a A fi V a A~<VpA 1 + \{A a A a + \). (1) 

Here, the Cj are dimensionless coefficients, and A is a 
Lagrange multiplier that enforces the constraint 

AM„ = -1. (2) 
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The Lagrangian ([I]) can be thought of as the low-energy 
description of a theory in which boost invariance is spon- 
taneously broken by the expectation value of A^, while 
spatial rotations and translations remain unbroken. The 
fixed-norm constraint eliminates the "radial" degree of 
freedom in field space, which is typically a ghost. We 
assume that is "minimally coupled" to gravity and to 
the rest of matter, so the total action is of the form, 



S = 



M 2 P 



J d i xy/^g[R + C A }+ J d A x^f^)C ri 



(3) 



Here Mp is the reduced Planck mass, and C m is the 
Lagrangian of ordinary matter, which we assume does 
not contain couplings to the aether field. 

The gravitational equations involve the 
energy-momentum tensor of the vector, 



T, 



ft V 



(-1/^)(SS A /Sg^). This is given by 



— V CT A v } - J^A^ - J^ v ^A"j + Y ttv + 
+ l -g llv C A + \A ll A v -c i A a AP{V a A lJi ){VpA v ), (4) 
where 

J a a = Cl V a A a + C2^V p A p + c 3 V a A a - c 4 A a A?VpA a , 

(5) 

and 

Y aP = a [(V 7 A Q )(V^ (3 ) - (VcAXV^)] . (6) 

Variation of the Lagrangian density |lj with respect to 
A leads to the field equation 

V a ( J a p) + c 4 A tt (V Q ^)(V^ 7 ) = \Ap, (7) 

whilst variation of the Lagrangian density with respect 
to the Lagrange multiplier A imposes the fixed norm con- 
straint ^j. 

The coefficients q are subject to both theoretical and 
phenomenological restrictions, which we collect in Ap- 
pendix [A] and summarize in Table [T] Their magnitude, 
relative to the symmetry breaking scale, can be esti- 
mated from dimensional analysis. The field redefini- 
tion A^ = Afi/M leads to the fixed norm constraint 
A^A^ = —M 2 , from which we may interpret M as the 
scale at which Lorentz symmetry is spontaneously bro- 
ken. In terms of the coefficients £j that would mul- 
tiply the action for the rescaled field A, the original 
coefficients are given by 01,2,3 = (M/Mp) 2 ci.2,3 an d 
c 4 = (M/Mp) 2 M 2 C4. We expect the dimensionless £1,2,3 
to be of order one, and the dimensionful £4 to be of order 
M~ 2 , which leads to 



M 2 
M 2 ' 



(8) 



For convenience, in what follows we use the abbreviations 



Cl3 = Ci+C 3 , C14 

a = ci + 3c 2 + c 3 , (3 



ci + c 4 , 

Cl + c 2 + c 3 . 



(9a) 
(9b) 



Note that a = 3/3 — 2ci3, so these abbreviations are not 
supposed to be an independent parametrization. Note 
also that our coefficients Cj and those of other works in 
the aether literature may have opposite signs. 



Cosmological dynamics 

Let us consider the dynamics of a spatially flat unper- 
turbed Friedman-Robertson- Walker universe in the pres- 
ence of the aether. Homogeneity and isotropy constrains 
the form of the metric and of the aether. With the line el- 
ement given by ds 2 = a 2 {rf) [— drj 2 + dx 2 ~\ we have, from 
Eq. 



A^ = (a" 1 , 0,0,0). 



(10) 



Substituting into the expression for the energy- 
momentum tensor Q, we find that the energy density 
and pressure of the vector field are respectively given by 



3a 2 
PA = 1C r , 2 rL , p A = - 



l6nGa 2 



{U 2 + 2U'), (11) 



where G = l/8irMp, % = a' /a and a prime denotes a 
derivative with respect to conformal time. Thus, Ein- 
stein's equations read 



H = 
W = 



8ttG 



3 

A-kG, 



cos 2 

— a p, 



cos 2 



(12a) 
(12b) 



where p and p are the energy density and pressure of the 
remaining matter fields (aether excluded) and 



G c . 



1 - 



G. 



(13) 



A comparison with the same equations in the absence 
of the aether shows that the effect of the vector field is 
merely to "renormalize" the value of Newton's gravita- 
tional constant [53] : the energy density and pressure of 
the vector field mimic that of the remaining components 
in the universe. 

On the other hand, the gravitational field created by 
isolated bodies is not exactly the same as that of Gen- 
eral Relativity, and in that sense the aether is a bona- 
fide modification of gravity. To lowest order in a post- 
Newtonian expansion, the potential sourced by a static 
and spherically symmetric body satisfies the Poisson 
equation Acf> = AttGnP, but with a modified gravita- 
tional constant [57] 



G 



N 



1 



Cl4\ 

2 ) 



G. 



(14) 



Hence, the aether also renormalizes the gravitational con- 
stant measured in "local" experiments, but by a different 
amount than in the cosmological case. Post-Newtonian 
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Condition 


Constraint 


Equation 


Solution of Einstein's equations 
Stability of Tensors 
Stability of Scalars 
Stability of Vectors 
PPN Limits 

Big-Bang Nucleosynthesis 


a < 2 
C13 > -1 
-2 < ci4 < 0, < 
2d < c? 3 (l + ci 3 ) 
see Equation 
ci4 + «<0.2 




16 

28 
38 
75 
Al 
A3 




Cherenkov radiation (assumes subluminality) 


see Equation 




A4 




Superluminal Tensors 
Superluminal Scalars 
Superluminal Vectors 


cis < 

(2 + ci 4 )/3< (2-a)(l + ci 3 )ci4 
2c 4 > -c? 3 /(l + Cl 3 ) 




27 
37 
74 




Anisotropic stress of long wavelength adiabatic modes 
Non-growing scalar isocurvature modes 
Subdominant contribution of vectors to CMB 


|ci 3 | < 1 

Q/ci4 > —1 

cY < c\ 


( 


58 
43 
118 





TABLE I: Summary of the theoretical and phenomenological conditions on the parameters of aether theories. We use the 
abbreviations a,/3,ci3 and C14, which are related to the standard aether parameters d through Eqs. |9| 



corrections lead to further deviations of General Relativ- 
ity, which place severe constraints on the aether param- 
eters. A summary of these and other constraints is given 
in Appendix [XJ Nucleosynthesis, in particular, places a 
bound on the relative magnitude of the two Newton con- 
stants, of the form [34] 



G 



- 1 



JY 



< 10%. 



(15) 



Note that, for positive matter energy density and pos- 
itive Newton's constant G, Eq. (12a) can only be solved 
if 3 



a < 2. 



(16) 



Remarkably, this condition does not follow from any of 
the perturbative stability arguments which we shall con- 
sider below, but merely from the existence of a cosmo- 
logical solution with positive energy density for ordinary 
matter. Note also that the Lagrange multiplier has a fi- 
nite value along the cosmological solutions. Contracting 
the vector field equations of motion (|7| with A@ we have 



It is well-known that this potential leads to power-law in- 
flation (35] . with a constant equation of state parameter 
w = Pip/Pip determined by the coefficient /i in the expo- 
nential. With a constant equation of state w the solution 
of Eqs. (|12a|) and (|12b| is then 



where 



with 



-H'/(aH 2 



1 



1 + 3w 



e-1' 



-/' 



(19) 



(20) 



is the conventional slow-roll parameter. Note that if 2 — a 
is sufficiently small, inflation may be de Sitter like even 
if [i is of order one. This broadens the class of "natural" 
inflationary models that do not require particularly flat 
potentials, though we shall not explore this possibility 
here. 



Cosmological Perturbations 



A 



{m 2 - C2W) 



(17) 



For later reference, let us consider the case where the 
matter sector consists of a scalar field with an exponential 
potential, 



C m = -^d^tpd^tp - V exp 



<P 



(18) 



3 We could have a > 2 if we allow G < 0. However, this leads to 
instabilities in the tensor modes, as we shall discuss in Section 
III. 



The background vector field ( 10 ) preserves rotational 



invariance, and so it is still convenient to use the stan- 
dard decomposition of perturbations in scalars, vector 
and tensors under spatial rotations: 

ds 2 = a 2 {ri) [-(14- 2(j))dr] 2 + 2{B * + S i )dridx i + 



+ (Sij - 2ip8 i:j + E tij + F hj + Fj.i 4- hij) dx i dx J 



(21) 



and 



-+SA , A^-ftCj + Vi-Si). (22) 
a a 
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Since the metric and vector fields are related to the La- 
grange multiplier by Eq. |7]), we also need to perturb the 
Lagrange multiplier, 



A = A + SX, 



(23) 



where Aq is the background value, given by Eq. (17) 



Variation of the second order action with respect to 5X 
leads to the linearized form of the constraint ([2| , 



5A U 



(24) 



the lifetime of the vacuum is then infinitely short, which 
makes the theory unviable. In a non-Lorentz invariant 
theory, the decay rate may be finite, and the vacuum 
may be sufficiently long-lived (see for instance [37]). In 
our case, Lorentz invariance is spontaneously broken, and 
the effective theory we are using is supposed to be valid 
only well below the symmetry breaking scale M. The 
decay rate is UV sensitive, so strictly speaking it is un- 
clear whether the theory can be made sense of in the 
presence of ghosts. However, to be conservative, we shall 
systematically exclude from parameter space the cases 
when ghosts are present. 

The primordial spectrum of tensor modes seeded dur- 
ing inflation is immediately obtained from (26), and is 



Here, <f>, B,ip, E,C are scalars, Si,Fi,Vi are transverse 
vectors, and hij is a transverse and traceless tensor. Note 

that hij and V, are gauge-invariant. Scalars, vectors and inversely proportional to their propagation speed, 
tensors decouple from each other in the linearized the- 
ory, so we consider each sector separately. In momentum 
space, our convention for the Fourier components is 



Vh(k) 



1 H 2 



f k ( V ) = f(r,,k) = 



d 3 x 



(2tt) 3 / 2 



/fa,x)exp(-ik.x). (25) 



III. TENSOR PERTURBATIONS 



7T 2 C t M 2 



(29) 



c 3 k=U 



Hence, the amplitude of the primordial tensor modes dif- 
fers from that in general relativity (for the same values 
of H and M P .) 



As discussed in [21] the presence of the aether modifies 
the propagator and the dispersion relation of the tensor 
modes. Substituting Eq. (21) into the action ([3]), with 
matter Lagrangian given by ( 18 ), expanding to quadratic 
order in hij and using the background equations of mo- 
tion we obtain, 



M 2 a 2 



[(1 + c 13 )h' • h' - d t h ■ 8 l h\ , (26) 



where h stands for a matrix with components hij and 
the dot indicates contraction of both indices (with the 
Euclidean metric). On short (subhorizon) scales, gravity 
waves propagate at a speed [35] 



4 



i 



(27) 



1 + Ci 3 

Classical stability of tensors thus imposes the condition 



1 + ci 3 > 0, 



(28) 



since, otherwise, high frequency modes grow exponen- 
tially fast. 

In the previous section we noted that the background 
solution only exists for a < 2, implicitly assuming that 
the "bare" Newton's constant is positive G > 0. Here, 
we note that for Mj, < 0, the coefficient in front of the 
kinetic term of hij has the "wrong" sign, and the two 
independent transverse and traceless tensor modes are 
ghosts. 

A theory with ghosts is quantum mechanically unsta- 
ble. The vacuum can decay by emitting positive energy 
particle plus negative energy quanta while conserving en- 
ergy. In a Lorentz-invariant theory, the phase space avail- 
able for the decay of the vacuum would be infinite, and 



(22 



aet 



IV. SCALAR PERTURBATIONS 

The scalar sector of Einstein-aether theories consists 
of the five scalars (f>, ip, B, E, C defined in Eqs. (21 ) and 
4 Thus, the aether enlarges the scalar sector by the 
ler perturbation C. 
It is convenient to introduce a gauge-invariant descrip- 
tion of the dynamical degrees of freedom. To this end, fol- 
lowing J57J [25] we note that the scalar part of the aether 
field can be represented by means of an auxiliary 
scalar field T through the identification 

A = ZLw 

M ~ {-r, v T?yi 2 ' 

where it is assumed that the gradient of T is everywhere 
time-like. Surfaces of constant T define a foliation of 
space-like surfaces, and we can think of T as a time 
variable. Since the background is aligned with the 
FRW temporal coordinate, the background field is given 
by T — T{r}). The perturbations (5T(t?,x) lead to the 
linearized spatial components Ai = —{a/T')57 From 
Eq. (22) we have Ai = adi(B + C), so it follows that 



ST 
V 



-(B + C). 



(30) 



In addition to the Einstein-aether sector, we must also 
include the matter sector. When the dominant matter 



4 The perturbation in the Lagrange multiplier <5A disappears from 
the Lagrangian after substituting the constraint to which it leads. 
To linearized order, this constraint is 5A° = —ij>/a, which we use 
to eliminate the scalar &A° in favour of the potential <f>. 
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component is the inflaton field ip, a convenient set of 
gauge-invariant variables is given by: 



C a = iJj-H(B + C), 
SN= '^5p + H{B + C). 



(31a) 
(31b) 



Geometrically, these can be interpreted as follows (see 
Fig. [Ij) Using (30 1, it is clear that the variable £ a is 



the curvature perturbation on surfaces of constant field 
T (i.e. on hypersurfaces orthogonal to the aether field 
A^). From the definition of Co. and SN it also follows 
that 



C = a + SN 



(32) 



is the curvature perturbation on surfaces of constant in- 
flaton ip. At the end of inflation and afterwards, C will 
describe the curvature perturbation on hypersurfaces co- 
moving with matter (excluding the aether.) On the other 
hand, 



5N = H 



Sip ST 
~ T 7 



HSt], 



(33) 



where 5r) is the amount of conformal time separating the 
surfaces of constant <p from the surfaces of constant T . 
Hence 5N can be interpreted as the differential e-folding 
number between these two types of surfaces. The velocity 
of aether with respect to the matter is given by 



(34) 



Hence, we can also think of the isocurvaturc perturbation 
W~ 1 6N as a velocity potential for the aether with respect 
to matter. 

In what follows, we consider the case of an exponen- 
tial inflaton potential, Eq. (18). This somewhat simpli- 



fies the analysis because the background solutions have 
a constant equation of state parameter p — wp. 

In addition, the behaviour of long wavelength pertur- 
bations of such a scalar field can mimic those of radia- 
tion and matter dominated eras for w = 1/3 and w = 
respectively. The "equivalence" applies only on large 
scales, because scalar perturbations and fluid perturba- 
tions have different sound speeds. Nonetheless, in Ap- 
pendix D we derive the form of the long wavelength adi- 
abatic and isocurvature scalar modes for generic matter 
content and expansion history. 



A. Short wavelength Lagrangian and stability. 

In Appendix[C]we discuss the Lagrangian for the scalar 
sector, and its reduction to a set of two gauge-invariant 
degrees of freedom (£ a ,5N) given by Eqs. (31a 31b). In 



the short wavelength limit, this Lagrangian is 



C 



2Z 



N 



[ { 5N)>z~e{5NY} + ^{Ca 2 -clkHl) + - 



(35) 



71 




<p = const. 



x 



FIG. 1: Geometrical interpretation of different perturbation 
variables. On hypersurfaces of constant inflaton ip, the cur- 
vature perturbation is while on hypersurfaces of constant 
aether T the perturbation of the spatial curvature is C, a . In 
the presence of isocurvature modes, both hypersurfaces do 
not agree. Their distance in conformal time is the variable 
5n, which measures departures from adiabaticity. 



where the ellipsis denotes terms which are subleading in 
the momentum expansion, and we have introduced 



<LrG c 



J N 



-2irc 2 t f3G c 



and 



G C os p 
Gn C14 



(36) 



(37) 



Here we have also introduced the slow roll parameter 
e = (3/2) (1 + w) and Gn as given in Eq. (14). For con- 



stant scale factor a, the residue Z a and sound speed c a 
agree with the corresponding quantities in a perturbed 
flat space, as discussed in [36], 5 Quantum stability re- 
quires Z a > 0, and classical stability requires c\ > 0. 
Recalling that stability of tensors demands c| > 0, and 
that ( 16 1 requires G cos > 0, we are led to the conditions 



2 < c 14 < and /3 < 0, 



(38) 



which in turn guarantees Gn > 0. (The case C14 = is 
singular, and has to be treated separately.) 

From Eq. ( 35 ) we can read off the normalization of 
the positive frequency modes associated with the "in" 



The above expressions are singular for = 0, but it is easy 
to show, following the derivation in Appendix [C] that f a is not 
dynamical in this case. 
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vacuum in the limit k\r)\ — > oo, corresponding to wave- 
lengths well within the horizon. The two independent 
mode functions are given by 



& ] -> 0, 



N e 



2k 



(a) 



Za 



1/2 —ic a krj 



a \ / 2c a k 



SN {a) -> 0. 



(39a) 



(39b) 



In the first mode, where ( a —> 0, the surfaces of constant 
aether field coincide initially with the so-called flat slic- 
ing, and SN is the number of e-folds separating the sur- 
faces of constant inflaton field from the flat slicing. This 
mode survives in the limit when there is no aether field 
(since one can still define the flat slicing surfaces). Hence, 
we may call this the inflaton perturbation. In the second 
mode, where SN —> 0, the inflaton is initially aligned 
with the aether, so that there is no inflaton perturbation 
in the aether frame. This mode survives in the flat space 
limit even when there is no inflaton field. Hence, we call 
this the aether perturbation. This can be thought of as 
one of the Goldstone bosons of the spontaneous Lorentz 
symmetry breaking. 

In the previous discussion we have assumed that the 
action ^ gives an accurate description of the aether 
up to sufficiently high momenta, so that ( 39 ) applies to 
scales well within the horizon. Hence, we require that 
the Einstein-aether as an effective theory should be valid 
at least up to some spatial cut-off A ^> H, where H is 
the Hubble rate during inflation. We expect the correc- 
tions introduced by the unknown physics above the cut- 
off scale to be at most of order of H/A (see for instance 

EE]). 

It should be noted that, on large scales, the fluctua- 
tions due to the aether will mix with those due to the 
inflaton. Hence, while in single field inflation perturba- 
tions have to be adiabatic, in Einstein-aether theories 
there should exist additional non-adiabatic modes, as we 
discuss next. 



B. Long wavelength modes 

The full Lagrangian in terms of the gauge-invariant 
variables £ a and SN is somewhat cumbersome away from 
the short wavelength limit, because the two modes are 
no longer decoupled. However, for long wavelengths the 
Lagrangian can be easily obtained from ( C8 ) and diago- 



nalized, but now in terms of a new pair of gauge- invariant 
variables (£, SN), where £ is the curvature perturba- 
tion on hypersurfaces of constant inflaton field, which 
we shall also refer to as comoving hypersurfaces, defined 
in Eq. (32). The long wavelength Lagrangian is given by 



C 



C c + C 



SN 



(40) 



Here the ellipsis denotes subdominant terms and the 
dominant ones are given by 



a 



(41) 



and 



n9 a 2 (kri) 2 



where we have introduced 



k = — 6 [ 1 



(SN)' 



1 



;(SNf 



cuj (l + 3w) 2 



(42) 



(43) 



As we shall see, there are a total of four independent 
long wavelength modes, which we derive in Appendix 
|B 2| Two of them have the property that SN — 0. For 
these, matter and aether are mutually at rest, and so we 
call these modes adiabatic. The other two have SN ^ 
and £ = 0, so we call them isocurvature, since there is no 
curvature perturbation on co-moving hypersurfaces. 



1. Adiabatic modes (5N = 0) 

In standard single field inflation, the non-decaying so- 
lution for the "adiabatic" perturbation £, which we de- 
note by Ci, stays constant on superhorizon scales. In Ap- 
pendix [D] we show that the same is true in the presence 
of the aether: 



£i = const. 



(44) 



for any expansion history (including the case where the 
equation of state changes abruptly in time). The cor- 
responding gravitational potentials in the longitudinal 
gauge are given by Eqs. (B19a): 



3(1 + w) 
(5 + 3w) 



4 Ci, 



Cl3 <?t Cl- 



(45a) 
(45b) 



The form of the two adiabatic modes (non-decaying and 
decaying) for an arbitrary expansion history and matter 
content is derived in Appendix|D] The decaying adiabatic 
mode is given by ( D8 ) , and it is characterized by 



4>2 = 1P2 oc Ha 
(2 = 0. 



(46a) 
(46b) 



It is worth mentioning that, although these adiabatic 
modes have the properties described in [39], they do not 
share the properties postulated in [40 42 . In particular, 
for the first adiabatic mode Ci, the anisotropic stress is 
non- vanishing (0i ^ ipi). 



2. Isocurvature modes (5N ^=0,(^ — 0) 



C. Power spectra 



As shown in Appendix |B 2[ for the case where the back- 
ground equation of state parameter w is constant, the two 
isocurvature modes behave as powers of conformal time: 



<5iVoc (-riY, 
where the exponents t are given by 

t± = 




(47) 



k, (48) 



and k is given in Eq. (43). Note that for k > 0, there 



is always a growing isocurvature mode. If we don't want 
this mode to grow out of control, then k should not be 
too large and positive, 



OO < K <Si 1. 



(49) 



In the following subsection we shall be more precise about 
the upper limit of this range (after discussing the over- 
all normalization of the corresponding power spectrum). 
Note that for a — — C14, we have n — and the dominant 
isocurvature mode stays constant on large scales, just like 
the adiabatic one. Hence, from the point of view of ob- 
servability of isocurvature modes, the interesting range 
of parameters is around a w — C14. 

From Eq. (B22c|, the gravitational potentials for the 



isocurvature modes are given in terms of SN by 



if} = -cia c 2 SN. 



and 



iP-<j> 



1. 



(50) 



(51) 



Hence, isocurvature modes have sizable anisotropic 
stress. It is also straightforward to check from the re- 
lations in the Appendix |B 2| that for the long wavelength 
isocurvature mode, the velocity of the aether with respect 
to matter is given by 



= n- 1 SN n = C 7 2 c, i 



(52) 



where in the first equality, we use Eq. (|34| and C is 
the scalar aether perturbation in the longitudinal gauge. 
From Eq. (Dll ) in AppendixlD] it is clear that at the time 



of a sudden transition in the equation of state parameter, 
the variable C and its derivative remain continuous for 
the long wavelength isocurvature mode. This means that 
the velocity field matches trivially: 



0. 



(53) 



where the square brackets indicate the discontinuity at 
the time of the transition. On the other hand, since 
the pressure changes abruptly at the transition, so does 
H' , and therefore the matching conditions for SN are 
[SN] = 0, [SN'] = (3/2) [w] nSN. 



As shown in Subsection |IV A[ the variables ( a and 
SN are uncorrelated on subhorizon scales. Hence, from 



Eqs. (32 1 and (39 1, it is clear that, at short wavelengths, 



the power spectra associated to SN and £ are given by 

2 



N 



{2-Kf 



z a + z_ 



N 



(54a) 
(54b) 



(2n) 2 

These spectra are valid for kr\ 3> l. 6 
1. Adiabatic modes 



For wavelengths comparable to the cosmological hori- 
zon, SN and £ are coupled to each other, and their evolu- 
tion will not have a simple form. Nonetheless, the evolu- 
tion of SN and £ is again simple in the long wavelength 
limit, as we saw in the previous subsection. In partic- 
ular, £ stays constant at long wavelengths. The power 
spectrum for £ will be approximately equal to its value 
at the time of horizon crossing, which we can estimate 
from (54a I by setting k/a = H. 



J N 



27T) 



H 2 



(55) 



where H 2 is evaluated at the time of horizon crossing. 
From Eqs. (36) and (37) we have 



J N 



Z a =[l 



Pec? 



J N • 



(56) 



Note that Z a is parametrically suppressed with re- 
spect to Zn by one power of aether parameters 
Cj ~ (M/Mp) 2 <C 1 and by one power of the slow roll pa- 
rameter e. Hence 



r c (k\ v \ « 1) 



3r 



[l + 0((3ec 2 )] 



(57) 



where p the energy density and ry^ is the time of hori- 
zon exit during inflation. Up to the small corrections 
introduced by the fluctuation of the aether, which are 
controlled by Z ai this expression is the same as the one 
in Einstein gravity, with Newton's constant G replaced 



6 Here, and for the rest of this section, we shall assume that the 
speed of propagation of aether is larger than or comparable to 
1. This is convenient so we do not have to introduce the scale of 
sound horizon crossing in the discussion of the adiabatic mode. 
Also, this assumption avoids the need of imposing the constraints 
due to Cherenkov radiation discussed in the Appendix. 
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with the effective Newton's constant G cos which appears 



in the Friedmann equation (|12a 

In summary, 
rameters a 



due to the smallness of the aether pa- 
the spectrum of primordial adiabatic modes 
does not change significantly with respect to the case of 
standard Einstein gravity. As we saw in the previous 
subsection, in the presence of the aether the adiabatic 
modes do not have the properties generally attributed to 
adiabatic perturbations. In particular, from Eq. ( 45b| , 
on super-horizon scales the non-decaying adiabatic mode 
has a non- vanishing anisotropic stress 



adiab 



^adiab 



Cl3 



^adiab 



(58) 



both in the matter and radiation dominated era. It is 
easy to see from Eq. (B6) that for a + cn — 0, the aether 



perturbation C behaves exactly like a massless field which 
propagates at the speed c a . Hence C oscillates while 
its amplitude decays as the inverse of the scale factor, 
C oc (l/a)e- iCakv . It then follows from (pi that (j> - tp 



also decays in inverse proportion to the scale factor, and 
hence it is suppressed by a factor of a(tk) / a(t Q ) , where 
tk is the time of horizon crossing. For modes which cross 
the horizon during the matter era, this means that the 
effect is suppressed with distance as 7 



ip-(f> 



C13 (fad) 2 {teq < k 1 < t , 



~ -C14), 

where we adopt the standard convention a(to) = 1. On 
the other hand, for modes that crossed the horizon dur- 
ing the radiation era, the scaling is with k^ 1 . This is in 
agreement with the result that at small scales the post- 
Newtonian parameter ip/(j> equals one, as in general rel- 
ativity 133]. Nonetheless, as a matter of principle, there 
could still be a distinct phenomenological signatures in 
the adiabatic sector imprinted on large scales. 

Constraints on the ratio ip/<f> on cosmological scales 
have been derived under several different assumptions, 
using combinations of different large scale structure 
probes [UHS]- At present, however, the constraints are 



7 Assuming that the aether parameters are small, these conclu- 
sions are easily extended to the case a 7^ — C14. In this case, 
Eq. §B6\ can be solved as the sum of the "homogeneous" equa- 
tion which is obtained by ignoring terms proportional to <f>, plus 
the contribution of a particular "inhomogeneous" solution. The 
first one takes the form Ch oc a~ (.l+d/2) e ic a kri ^ w h ere 



- c 13 c? 1 + 



(59) 



This leads to 

ij>~<j> 



~ C13 [(fct r( 2+d > + o(fc*or 2 ] iu q < k' 1 « t ). 

(60) 

This applies to modes that entered the horizon during the matter 
era. For those which crossed the horizon during the radiation era, 
the scaling is with one less power of k in the denominator. 



quite weak, and it appears that values of ip/4> of order 
one are still consistent with the data. 



2. Isocurvature modes 

Next, let us consider the spectrum of long wavelength 
isocurvature modes Vsn- The phenomenological situa- 
tion depends on whether a + C14 is positive or negative. 
If a < — C14, then SN decays on superhorizon scales, 
during and after inflation. Hence, these modes will re- 
main insignificant with respect to the adiabatic ones. If 
q = — C14, then there is a constant non-decaying isocurva- 
ture mode, and SN stays constant on superhorizon scales. 
Finally, for a > —cn there is a growing mode and SN 
can be very large at the time of re-entry even if it was 
small at the time of horizon exit. 

Phenomenologically, the most interesting case seems to 
be the limit |a + ci4| <C | C14 1 , in which the supercurvature 
mode SN stays approximately constant on large scales. 
Otherwise, either the mode is too suppressed to be of 
any significance, or it grows too fast to be compatible 
with observations. In this case, the exponent t± for the 
dominant mode can be approximated by 



t ■ 



l + 3w 
5 + 3w 



(62) 



where k is given in (|43|, and we have \i\ <C 1 both during 
inflation and afterwards. At the time of horizon crossing, 
the adiabatic and isocurvature modes have comparable 
amplitudes, V ~ (2ir)~ 2 Z^H 2 , and these will remain 
roughly comparable throughout cosmic history up to the 
present time provided that t is sufficiently close to zero. 
In order to assess how small it would have to be, we can 
make a rough estimate of the evolution of the amplitude 
of SN from the time of horizon crossing during inflation 
to the time of equality: 



UN 



Vcq 
?7rh 



z]i 2 He^~^ N < 1. 

(63) 

Here, the subindices i and r refer to inflation and ra- 
diation era respectively. Assuming Z]^ 2 H ~ 10~ 5 , as 
follows from the normalization of the adiabatic modes, 
we find that for 



i r - it fa 1 n r < ln(10 5 )/iV 



(64) 



the perturbation SN remains within the linear regime up 
to the time of equality of matter and radiation. Here, we 
have neglected Ki , which is suppressed with respect to n r 
by a slow roll factor (we are assuming that the aether 
parameters are the same today than they are during in- 
flation), and 



N ~ 60 



(65) 



is the number of e-foldings of inflation after the mode 
with co-moving wavenumber k ~ rj cq first crossed the 
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horizon. Note also that, according to (50 1, the contribu- 



tion of the isocurvature mode to the gravitational poten- 
tial is suppressed by C13, 

Vw = -cis <? t < 10- 5 , 

where the last inequality is the observational bound on 
the gravitational potentials. For C13 c 2 < 10 -5 , V'isoc can 



remain small enough even if the inequality (64) is satu- 
rated, so that SN ~ 1. Also V'isoc becomes comparable 
to the contribution of the adiabatic mode i/Viiabj when 
the inequality K < — 31n|ci 3 c 2 \/AN is saturated, and 
combining with ( 64 ) , we require 



< — min{ln(10 5 ), 



•ln|ci3 c 2 |}. 



(66) 



Let us estimate what the physical implications of the 
isocurvature perturbations might be. On one hand, they 
would induce maximal anisotropic stress on large scales, 
as can be seen from Eq. (51 ), 



ipii 



1, 



which means that there would be a sizable difference 
between the two gravitational potentials <f> and pro- 
vided that the contribution of the isocurvature mode is 
comparable to that of the adiabatic mode. For n = 
we have Vw ~ c 13 cjSN ~ Ci 3 cf<0adiab 5 and so from 



(58) both adiabatic and isocurvature modes contribute 



to the anisotropic stress in a similar amount (unless c 2 
is very large). However, if k is small and positive, then 
the isocurvature mode grows on superhorizon scales, and 
will contribute more to the anisotropic stress than the 
adiabatic one. As argued in the previous subsection, the 
difference 1 — (</>/V0 decays after horizon crossing, and 
so does its magnitude as a function of co-moving scale, 
which roughly goes as k~ 2 for modes which crossed the 
horizon during the matter era, and as fc _1 for modes that 
crossed before the time of equality [see Eqs. (60), (61 )]. 



Another possible signature might be due to preferred- 
frame effects due to the motion of matter with respect 
to the aether [48j E2] , such as a dipole anisotropy in the 
gravitational potential of massive bodies. The primor- 
dial perturbations cause the aether to point in different 
directions at different places in the observable universe. 
Hence, the velocity of matter with respect to the aether 
(and the corresponding gravitational dipole, for instance) 
would have a random distribution. From ( |34[ ) an isocur- 
vature perturbation with wave-number k, induces a rel- 
ative speed of the aether with respect to matter given 
by 

v = -SN. 

n 

When the mode reenters the horizon during the radiation 
or matter era, at time ~ a/k, we have 

v~6N(t k ). 



This has to be compared to the peculiar velocities in 
bound objects at the same scale, which is of order 
C, 1 / 2 ~ 10~ 3 . Hence, the effect of the peculiar velocity 
of the aether will be subdominant unless SN has grown 
from the time of horizon exit, in such a way that at the 
time of reentry it is at least of the order C 1 / 2 . This pos- 
sibility exists, since we have seen that SN has a grow- 
ing mode for a > — C14. Because of that, the velocity 
of the aether at the time of horizon crossing could even 
approach moderately relativistic speeds without compro- 
mising the validity of the linear approximation and with- 
out contradicting current observations. [Note from Eqs. 
(50) and (51 ), that the gravitational potentials along the 



isocurvature mode, and hence their effects on the CMB, 
are suppressed by a factor C13, which can be very small]. 

The velocity field of the aether is strongly correlated 
with the amplitude of adiabatic modes, since both have a 
common origin in the amplitude of the short wavelength 
mode SN when it first crosses the horizon during infla- 
tion. Should the velocity field of the aether be detected, 
such correlation would indicate that the velocity field has 
a primordial inflationary origin. 

We may define a power spectrum V v for the longitudi- 
nal velocity field of the aether through the equation 

9ir 2 kk 
(v^k^X)) = -p-P o fa,fc)^5(k-]0. (67) 

Note that at the time of horizon exit during inflation, we 
have 



Vv ~ V 5N ~ Vr 



10" 



(68) 



where the last estimate follows from observations. How- 
ever, since the perturbation SN grows on large scales for 
< k <C 1, we can have V v ~ Vsn ^ V$ at the time 
of horizon reentry. As we shall see in the next section, 
vector perturbations can give an additional contribution 
to the velocity field (which can of course be disentan- 
gled from the scalar isocurvature contribution from the 
fact that the corresponding velocity field is transverse). 
It turns out that the scalar component and transverse 
vector component of the velocity field obey the same 
equation of motion on large scales. Hence, we defer the 
discussion of the spectral properties of V v on currently 
observable scales to the next section. 



V. VECTOR PERTURBATIONS 

In a universe dominated by a scalar field there are no 
vector perturbations. Perfect fluids do support vector 
perturbations, but they decay as the universe expands. 
By contrast, the aether contains a massless vector field 
(under spatial rotations), which renders vector metric 
perturbations dynamical. For a certain range of param- 
eters, these modes can grow on large scales, leading to 
potentially interesting signals, or to constraints on the 
parameters of Einstein-aether theories. 
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A. Short wavelength stability 

As in the case of the tensor modes, we can read off 
from the action for the vector perturbations whether the 
vector sector in Einstein-aether theories is both quantum 
and classically stable on short scales. Inserting the ex- 
pansions (21 1 and (22 1 into the action ([3]), expanding to 
quadratic order in the vectors S and V, and using the 
background equations of motion, we obtain the following 
Lagrangian for the vector perturbations 



L (2) 



M 2 a 2 



- C i 4 V' 2 + ^(l + Cl3 )a t Q-a i Q+ 

+ cidiV ■ d l V + ci 3 a 4 Q • d l V+ 

+ a(H 2 - W) V 2 + c 14 {H 2 + n')V 2 



(69) 



where we have introduced the gauge-invariant combina- 
tion 



Q = F' - S 



(70) 



(the vector perturbation V is also gauge-invariant). Note 
that Q is not a bona-fide Lagrangian variable, since its 
definition (70) relates it to the time derivative of F. 



Hence, we shall merely think of it as shorthand for the 
right hand side of ( 70 1 . Variation of Eq. ( 69 ) with re- 



spect to S gives the response of the metric to a given 
perturbation of the aether field, 



Q = -Cl3 



c 2 V. 



(71) 



In the canonical (first order) formalism, this equation 
corresponds to the vanishing of the canonical momentum 
conjugate to F, ITf = 0. Upon substitution of this con- 
straint back into the first order Lagrangian, one is left 
with a Lagrangian for the single canonical pair formed 
by V and its conjugate momentum Ilv. Rewriting this 
reduced Lagrangian back in second order form gives 



f(2) 



M 2 



ci 4 £' 2 + a(H 2 - H')i 2 



•ci 1 



r 2 r 2 

2ci 



(72) 



where, for convenience, we have introduced the rescaled 
variable 



& = aVi. 



(73) 



As in the tensor case, the absence of ghosts requires 
the coefficient in front of £' 2 in the Lagrangian ( 72 ) to be 



positive, and classical stability requires that the squared 
speed [36] 



C14 



2 Cl 



(74) 



be non-negative. Therefore, stability in the vector sector 
demands both 



C14 < and ci < 



-13 



2(l + c 13 )- 



(75) 



In a Minkowski background, the two modes in the vec- 
tor sector are massless fields, which we may interpret as 
two of the Goldstone modes of the broken boost invari- 
ance. The broken generators transform as a spatial vector 
under the unbroken group of spatial rotations, so the cor- 
responding Goldstone bosons transform as a vector. This 
can be decomposed into a transverse part and a longitudi- 
nal part. The longitudinal Goldstone (with helicity zero) 
is of course part of the scalar sector, which we discussed 
in the previous section. It should be noted that Lorentz 
invariance is generically broken in any curved spacetimc. 
For instance, if the spacetime curvature is non-constant, 
the gradient V M i? defines a non-zero vector field which 
is not invariant under Lorentz-transformations. What is 
particular about the Einstein-aether is that the breaking 
of Lorentz-invariance has physical consequences, namely, 
the existence of Nambu- Goldstone bosons, whose disper- 
sion relations approach non-relativistic expressions in the 
high-momentum limit, and whose masses vanish in flat 
space. 8 



B. Solutions during Power-Law Inflation 

Variation of ( [72| with respect to £ leads to the equation 
of motion for the vector perturbations, 

e:' + c ?,fc 2 6 + — {n 2 -n')i l = Q. (76) 

C14 

In terms of the original variable Vj, = £i/a, we have 



V" + 2UV( + c 2 k 2 V t + 
1 + ^W 

Cl4/ 



C i4 



Vi = 0. (77) 



In a universe that undergoes power-law inflation ( 19 1 this 
can be solved in terms of Bessel functions and we have 

2Mp y — C14 a 



Here, 



t+ - t- 




(79) 



8 In an arbitrary spacetime, the "mass" of these bosons is non-zero, 
as illustrated in the case of a FRW universe by the cont ribu tion 
to the effective mass of the last term in the Lagrangian |72|l. 
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where t± is given in ( 48 ) . The parameter q is defined in 
Eq. ( 19 ) and is a normalized transverse polarization 
vector, e • k = and e 2 = 1. For a given wave number k, 
there are two such linearly independent polarizations, or- 
thogonal to k. We have chosen the amplitude of Vi in Eq. 
( 78 ) so that the solution has the appropriate normaliza- 
tion of a positive frequency mode in the limit 77 — > — 00. 
Note that the factor q(q + 1) is non negative if the null 
energy condition is satisfied (w > —1). 

The long wavelength power spectrum of vector pertur- 
bations created during inflation is defined by 

(Vi{rt, k^k')) = ^TMM)n^(k-k'), (80) 

where = Sij — kikj/k 2 projects onto the subspace 
orthogonal to k. At the time of reheating ?y r h, we have 



n v = 3 — 2v, 



(81a) 
(81b) 



M 2 \c 14 q 2 J (2?r) 



tt) 3 V 2 J 



exp 



(81c) 



In Eqs. (81 ) is the value of the Hubble constant at 
the end of inflation, and k^ is the mode that crossed 
the cosmic horizon N e-folds before the end of inflation 
(Ifcjvwl = !)■ F° r the mode that is entering the hori- 
zon today, the value of N depends logarithmically on the 
unknown reheating temperature, and typically equals 50 
to 70 e-folds (see for instance [50].) It is important to 
realize that the time at which the spectrum is evaluated 
matters, since the vector modes do not freeze out at hori- 
zon crossing. The superscript "rh" is meant to imply that 
the power spectrum describes the amplitude of the modes 
just before the end of inflation. Likewise, we may define 
the spectrum of the corresponding metric perturbation, 
which according to Eq. ( 71 ) is given by 



a2 _ 2 2 a2 
■^Q — c 13 c e /l V 



(82a) 
(82b) 



In the limit in which de Sitter inflation is approached 
(q — > —1) the index v tends to 1/2, so the ampli- 
tude of long-wavelength perturbations is proportional 
to exp(— 2N). Hence, velocity perturbations on obser- 
vationally accessible scales are very small in this limit 
[29j |30] . On the other hand, in typical inflationary mod- 
els q differs from —1 and the rate of decay can be smaller. 
In fact, if 



a > 



-2c 



14 



-2ci4 



(83) 



the combination 3 — 2^ would be negative, and long- 
wavelength perturbations would be amplified exponen- 
tially with N, in stark contrast with the de Sitter case. 



It turns out, however, that we do not need to deviate 
much from v = 1/2 in order to have an observable sig- 
nal. As we shall see, even if the long wavelength velocity 
field is very tiny at the end of inflation, it may resurface 
from obscurity during the radiation and matter era, so 
that it can be quite sizable at the moment of horizon 
reentry. 

Indeed, the behaviour of long wavelength vector per- 
turbations is completely analogous to that of the scalar 
component of the velocity field Vi which we studied in the 
previous section. To see this, we note that in the long 
wavelength limit, Eq. (77) is the same as Eq. (DID for 



the isocurvature perturbation. The latter is written in 
terms of the variable C in the longitudinal gauge, which 
according to Eq. ( 52 ) is proportional to the longitudinal 
velocity field of the aether Uj. Hence, on superhorizon 
scales, longitudinal and transverse velocity fields satisfy 
the same equation of motion Eq. ( 77 1 . In particular 



Foe vex (k/a)n~ 1 6N oc 77 



l-q+t± 



(fcr?«l), (84) 



where t± is given in (48). For a — — C14, the velocity field 



decays exponentially during inflation. Nonetheless, as we 
saw in the previous section, the isocurvature perturbation 
SN stays frozen on superhorizon scales (except at the 
transitions where the equation of state changes, where 
the dominant mode changes also by factors of order one) . 
This can lead to a sizable velocity field v of order SN at 
horizon reentry. The overall normalization of v and V is 
different, but it is clear that the relative size of V and v 
at horizon reentry is determined by their relative size at 
the time when they exit the horizon during inflation. In 
other words, the spectra of long wavelength modes are 
related by 



Vy(y,k) = V v (r]k,k) 



C14 
Mj 



C14 



(85) 



where the relative normalization can be read off from the 
corresponding short wavelength actions, and e is the slow 
roll parameter during inflation. Note that e is of order 
of a few percent, while the Einstein aether parameters 
such as C14 ~ (M/Mp) 4 are suppressed by the square 
of the symmetry breaking scale over the Planck scale. 
Unless M and Mp are very close, we expect C14 <C e. 
Therefore, parametrically, we expect that the transverse 
vectors may give a much bigger contribution than the 
scalar isocurvature modes. For that reason, it is very 
important to assess their impact on observables such as 
the CMB, as we do in Subsection |VIB| 

One may worry that if the primordial amplitude of vec- 
tors due to quantum fluctuations generated at horizon 
crossing decays during inflation, then it may be insignifi- 
cant compared with the contribution of non-linear effects 
which source the vector modes at later times. However, 
in order to construct a vector from a quadratic expres- 
sion involving scalars and tensors, it is necessary to use 
at least one derivative. Because of that, the terms which 
may source the vectors from the scalar and tensor sector 
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are momentum suppressed, and hence they also decay in 
inverse proportion to the scale factor. We conclude that 
if the initial amplitude of the vectors at horizon crossing 
is sizable, compared to that of scalars and tensors, then 
we can safely use linear evolution in order to determine 
its amplitude at the end of inflation (even if that ampli- 
tude is exceedingly small). Vector modes can still grow 
during radiation and matter domination from that ini- 
tial tiny amplitude, so that their effect on cosmological 
observables may be important. 



VI. VECTOR CONTRIBUTION TO CMB 
TEMPERATURE ANISOTROPIES 

As we discussed above, the power spectrum of the 
transverse velocity field may easily dominate over that 
of the longitudinal component. It is therefore of inter- 
est to determine the imprint that this power spectrum 
may have on CMB observations, to which this section is 
devoted. 



Solutions During Radiation and Matter 
Domination 



In order to find the power spectrum of the vector 
modes at the time of recombination, we must first evolve 
it from the time of thermalization through the radiation 
and matter dominated epochs. For a set of perfect fluids 
which do not interact with each other, the conservation 
equation V^T/ 1 — holds for each fluid component, 
which we label by (fc). This leads to a homogeneous equa- 

( k) 

tion for the gauge-invariant velocity perturbation 5u\ 
that does not contain metric perturbations, 



a 3 {Pk +Pk)Su 



O) 



0. 



(86) 



from where it follows that Sui/a oc l/[a 4 (p+p)}. Eq. (86) 
tells us that if Siii = initially, then it will not be gen- 
crated as long as the perfect fluid description is valid. 
Furthermore, Sui/a decreases during cosmological evo- 
lution, except in the radiation dominated stage, where 
it stays constant. Hence, in what follows, we assume 



that Su 



for matter and radiation. Then, the only 



contribution to the metric perturbations stems from the 
aether, and it can be shown that the equations of motion 
in the vector sector are still given by Eqs. ( 71 ) and ( 76 1. 
The general solution of Eq. (76) during a stage of 



cosmic expansion in which a cx r\ q is proportional to a 
linear combination of Bessel functions. For our present 
purposes, it will suffice to work with the long-wavelength 
approximations 



A. 



(87) 



where v is given by Eq. (79) and ^ = £ • ej. For real 
values of v, A+ is the amplitude of the dominant mode, 
and A- is the amplitude of the subdominant mode. 

In order to determine the mode amplitudes A\ and A r _ 
during radiation domination, we simply demand that £ 
and its time derivative be continuous at a sudden tran- 
sition from inflation to radiation domination. We expect 
this approximation to be valid for scales much longer 
than the duration of reheating. Proceeding in this man- 
ner, and dropping the contribution from the subdominant 
mode we find that the amplitude of the dominant mode 
changes during reheating by a factor 

A\ v l + v r 



A\ 



2v r 



where the superscripts label the expansion epoch (i for 
inflation and r for radiation domination), and the sub- 
scripts label the different modes (+ for the dominant 
mode and — for the subdominant one.) 

Eq. (76) has an exact solution at long- wavelengths 



during radiation and matter domination, which we can 
use to determine the change in the mode amplitudes dur- 
ing the transition from radiation to matter domination. 
Since this change is typically of order one, we shall ne- 
glect it, and assume that the amplitude of the growing 
mode at the transition remains unchanged. 

Once a mode enters the "sound horizon", c v krj = 1, 
the field starts oscillating. In the limit c v krj 1, the 
solution of Eq. ( 76 ) that approaches the growing mode 
at early times is 

£(77) = A+ C osc cos [c v kri — ip] , where (89a) 
c 13 (v + l) (2 



Cy kTJ^. 



(89b) 



and ip is fc-independent phase. Note that the amplitude 
of the oscillations A+C osc is roughly the value of £ at 
horizon entry. 



Collecting then the results from Eqs. (73), (87) and 



( 88 1 we find that during matter domination the trans- 



fer function for the vector perturbations 1\, which we 
implicitly define by the relation Q(r?) = 7fc(?7)Q(f7rh) is 



Tk~Ty.l 



where 



^cq 
Orh 

a rh 

Cosc cos(c v krj) 



— , ki] cq < kr/ < c \ 

C™ c cos(c„fc?7), ki] cq < c^ 1 < kr), 

c" 1 < ki] oq < krj, 
(90a) 



T = 



2v r 



(90b) 



and in C% sc the transition time 77* equals Tjrh for p = r 
and ?7 q for p — m. The first line in Eq. (90a) holds for 
those modes that have not entered the sound horizon at 
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time 77. The second applies to those which enter between 
the time 77 oq of equality between matter and radiation 
densities and time 77, and the third one to the ones which 
enter between reheating and the time of equality. The 
power spectrum of Q at any time after reheating is given 
by 



■Pq(v) = \%(v)\ 2 V, 



Hi 



(91) 



B. Impact on the CMB 

We derive in Appendix [E] the contribution of vector 
perturbations to the temperature anisotropies in the cos- 
mic microwave background radiation. In order to deter- 
mine the angular power spectrum and relate it to the pri- 
mordial spectrum, it is convenient to rewrite Eq. (E21) 
in Fourier space, 



ToJo 



l'Q(?7dcc, k ) ex P (* k • l?7dcc) + 



»/0 

J dr}l-Q'(ri,k)exp(ik-lri) 



(92) 



The contribution of the two terms on the right-hand-side 
of (92) is similar to that of scalar perturbations. The 



first term is the analogue of the Sachs- Wolfe effect, which 
relates the temperature anisotropies to the state of the 
perturbations at last scattering. The second term is the 
analogue of the integrated Sachs- Wolfe effect, which takes 
into account the change of the metric potentials along the 
line of sight, and vanishes if the latter are constant. 

The angular power spectrum C/ is defined by the rela- 
tion 



'ST ST 
, Jo Jo 



2£- 



4tt 



-Pt(n-m), (93) 



In these equations x = k(rjo — 77), the jg are the spherical 
Bessel functions of the first kind, the primordial spec- 
trum is given by Eqs. (82) and the transfer function Tu 
by Eqs. (90). After an integration by parts, Eqs. (94) 
agree with the expression derived in [51] by somewhat 
different methods. Note that x is the ratio of the co- 
moving distance to time 7/ divided by the wavelength of 
the perturbation 1/fc. Thus, x is the inverse of the angle 
that an object of comoving size 1/fc at (comoving) dis- 
tance 770 — 77 would subtend on the sky at time 770. It will 
be useful to consider separately those modes that enter 
well before and well after decoupling. These contribute 
to the temperature anisotropies, respectively, on small 
and large angular scales. 



The structure of expressions ( 94 ) still reflects the two 



contributions to the temperature anisotropies we men- 
tioned above. The value of Af captures the analogue of 
the Sachs- Wolfe effect, while the value of the integrated 
term I captures the analogue of the integrated Sachs- 
Wolfe effect. For scalar perturbations, the Sachs- Wolfe 
effect is dominant except on the largest scales, because 
the gravitational potential remains constant until rela- 
tively recently. For vector perturbations however, this is 
not always the case. To see this, it is useful to realize 
that we can employ the same approximations developed 
to study the contribution of tensor modes to the tem- 
perature anisotropies |52j . We begin by noting that, for 
£ ^> 1, the Bessel function can be approximated by [S3] 



0, 



j e {x) « I cos y/y - (£ + 1) arccos (^~) 



x <£ + 



x > £ + 



x l/2 y l/4 

(95) 

where y = x 2 — (£ + 1/2) 2 . Because the integrand is 
negligible for x < £ + 1/2, only modes that have entered 
the horizon by today, Xd oc ~ krjo > £+1/2 can contribute 
to the temperature anisotropies on angular scales £ ^> 1. 



where n and rh are two directions on the sky, and the 
Pi are Legendre polynomials. Because scalar, vector 
and tensor perturbations are uncorrelated, their contri- 
butions to the temperature anisotropies ad d in quadra- 
ture, C e = C| + Cj + C\. Inserting Eq. Q into the 
left-hand-side of Eq. ( 93 ) , using Eq. ( 80 ) , and compar- 



ing to the right-hand-side of Eq. ( 93 ) we find after some 



work that the contribution of vector perturbations to the 
angular power spectrum is given by 



Af 



^£(£+1) J^V$\Af~l\ 2 , (94a) 

je(Xdec) 



Tk (r/dec) 



•^dc 



X = 



dx 



dTk jejx) 
dx x 



(94b) 
(94c) 



Large Angular Scales 

Large angular scales correspond to modes that cross 
the sound horizon after decoupling. Let us estimate the 
contribution of the integrated term T on those scales first. 
From Eq. (90a), the derivative of the transfer function 
dTk/dx is dominant during the interval Xk ^5 x < Xdcq, 
where Xk = k(rjo — r/k) corresponds to the time of sound 
horizon crossing rjk = l/(c v k). After horizon crossing, 
this function oscillates with period 2it/c v and a slowly 
varying amplitude. On the other hand, the ratio j(x)/x 
changes slowly in the interval Ax < 1 ^ I < x. As- 
suming that c v is not much smaller than 1, we have 
Ax = x^cq — Xk ~ kr/k — l/c v < 1, and we can pull the 
factor j(x)/x out of the integral. What remains is a 
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boundary term that can be readily evaluated, 



, dT k ji(x) j ( {x dcc ) 
dx — k, [7fc(r; d cc) - Tk{Tl*)\ 



dx 



(96) 

where the effective lower limit of integration 77* is of order 
Vk- 



The dominant term in the right hand side of Eq. ( 96 ) 



depends on whether the long wavelength modes are 
decaying (which happens for v m < 3/2), or growing 
(y m > 3/2) before horizon crossing. In the first case we 
have 



T ~ ji(Xdcc) rj- I \ 

-L « lk(Vdec), 

Xdcc 

while in the second case we have instead 

je(Xdec). 



1 



(97) 



(98) 



Therefore, comparison of Eqs. (97) and (98) with (94b) 



shows that AT and I are of the same order if there is no 
growing mode during matter domination (5 + v m < 2), 
and that X > M otherwise (y m > 3/2). 

We are ready now to calculate the angular power spec- 
trum on large angular scales, which are dominated by 
modes that crossed the vector sound horizon after de- 
coupling {c v krj dcc < 1). The relevant expression fo r the 
transfer function is given by the first line of Eq. (90a) 
at rj = r/dcc- If there is no growing mode during matter 
domination, the contributions from the integrated and 
non-integrated terms in (|94[ ) are roughly equal, and sub- 
stituting Eqs. (82) and (1911) into (94) we get 



Cj ^^A 2 Q T k 2 {iid CC ) J 



dx £{£ + 1) 



J'O 



x^jUxo), (99) 



where we have chosen fc/v in Eq. (81) to be the mode 
that is crossing the horizon today, kj^r/Q = 1, and used 
that Xdcc ~ x ee kr/o- 

From reheating to the time of decoupling, the ampli- 
tude of the vector modes changes by 7fc(r/dcc)- The spec- 
trum is thus proportional to the square of the transfer 
function times the primordial amplitude Aq . This factor 
is independent of angular scale, since we are taking the 
long wavelength limit. The angular dependence in the 
last equation can be estimated as follows. The Bessel 
function is negligible for xq < £, and rapidly decays at 
xq > £, so the anisotropics are dominated by xq ~ £. In 
the integrand, the maximum of the Bessel function is of 
order 1/xq, and the two remaining factors of xq in the 
denominator "cancel" the enhancement proportional to 
£(£ + 1) one would otherwise have. In summary, we have 
Cjoz£^- 2 . 

More precisely, since the period of oscillations of the 
Bessel function is much shorter that any other character- 
istic scale in the integrand of Eq. ( 99 1 , we may replace 
the oscillations with their average, 1/2. If the spectral 



index n v is not too blue (n v < 4) , the dominant contribu- 
tion to the integral is given by the value of the integrand 
at x s» £, so the angular power spectrum becomes (for 
v m < 3/2) 



£{£+l)C\ 



2nAf 



V 



2v r 



<^dcc 



2v T -\ 



cq 



u m -2u T - 



£ n -». 

(100) 



This expression is valid for those scales that entered the 
vector horizon after recombination, which corresponds to 
£ < 50/c„ (for a ACDM model with Q A = 0.7.) 

If there is a growing mode during matter domination, 
the integrated term ( 94c ) yields the dominant contribu- 



tion to the temperature anisotropies. Proceeding along 
the same lines as above, we find that in this case the 
angular power spectrum is (for v m > 3/2) 



£(£ + l)CY 



2-kA% 



2v r 



u, cq 



21A-1 



7 riu+3— 2i/ 



(101) 



Of course, in the crossover case v m = 3/2 the two angular 



power spectra ( 100 ) and ( 101 ) agree. 



Small Angular Scales 



For the scales that enter the vector horizon before de- 



coupling a precise estimate of the integrated term in ( 94c 
becomes more difficult. For these modes the derivative 
of the transfer function is an oscillatory function, whose 
amplitude decreases in time. Hence, it is most impor- 
tant at earlier times x sa x dcc and sharply decays within 
an interval Ax — fc^dec- Whereas the latter is small for 
modes that cross after decoupling, for those scales that 
enter the horizon well before that time Ax = kr)d cc is 
large, and the approximation of a constant Bessel func- 
tion that led to (97) breaks down. Nonetheless, if we are 



interested in the order of magnitude of the Bessel func- 
tion, and not in the oscillations, we can still use Eq. ( |97| , 
since the "amplitude" of je(x) only changes significantly 
within Air = Xdcc ~ kr/o ^ kr/d cc . In that case, the 
integrated term is at most of the same order of magni- 
tude as the non-integrated one, and Eqs. (94) imply that 



the temperature anisotropies at any given angular scale 
will depend on the vector anisotropies on the appropriate 
comoving distance at the time of decoupling. 

Under the assumption that Af and I are of the same or- 
der, the angular power spectrum on small angular scales 
can be now calculated as before. For simplicity, let us 
concentrate on relatively small scales, which cross the 
sound horizon before equality of matter and radiation 
densities. The relevant expression for the transfer func- 
tion is given by the third line of Eq. ( 90a ) . Substituting 



16 



Eqs. (j82| and ( |90a[ ) into Q we obtain 



d 4^^T^ 1 (r7dcc)x 
dxo + 1) 



-2v T - 



'cos 2 



Xq 



Vo 



xo jt(xo). 



(102) 



As before, the power is proportional to the primordial 
contribution x nv times an additional factor x~ 2l/ _1 , 
which just reflects that modes enter the horizon at dif- 
ferent times, and thus evolve differently. The cosine rep- 
resents a snapshot of the "acoustic oscillations" of the 
vector perturbations at decoupling. 

Before we proceed, we should mention an additional ef- 
fect that influences the anisotropics on very small scales. 
So far, we have been assuming that the decoupling of 
the photons from the baryons is instantaneous. This is 
an accurate approximation for scales in which the argu- 
ment of the cosine in Eq. 



(90a 



does not change much 
during the duration of decoupling. On scales in which 
the cosine does change significantly, the spread in time 
at which a photon last scatters dampens the fluctuations 
by an exponential factor exp (— Xq/2<t 2 ) 52J . A similar 
suppression is also due to Silk-damping, which originates 
from the breakdown of the tight-coupling approximation 
at scales of the order of the mean free path of photons 
in the plasma. For the observed values of the cosmologi- 
cal parameters, both effects yield an overall value of the 
suppression scale a « 500. 

Due to the exponential damping, the integral over 
modes converges for any power-law spectrum. As before, 
if the effective spectral index is not too blue, the domi- 
nant contribution to the integral is given by the value of 
the integrand at xq rj £, so the angular power spectrum 
becomes 
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(103) 



This equation is qualitatively valid at small angular 
scales, those corresponding to modes that crossed the 
horizon before equality, t > 120/c„. The acoustic oscilla- 
tions subtend an angle Cv^dec/vo on the sky, the ratio of 
the comoving size of the sound horizon at decoupling to 
the comoving distance to the last scattering surface. A 
plot of the angular power spectrum for vector modes for 
a specific set of parameters is shown in Figure [2j 



t((+l)C, 




50 100 



500 1000 



FIG. 2: The contribution of vector modes to the temperature 
anisotropy power spectrum for models with C14 = —a and 
c v = 1 (the normalization is arbitrary.) The black (contin- 
uous) curve is the added contribution of the integrated and 
non-integrated terms, Eq. ( 94 1 . The contribution of the in- 
tegrated term alone is shown in red (dashed-dotted), while 
the contribution of the non-integrated term alone is shown in 
blue (dashed). On large angular scales, the spectrum is well 
approximated by Eq. (101 1. On small scales, Eq. (103 1 gives 



a qualitatively correct approximation. 



that of the tensor modes [52j [54] , which, in the limit of 
instantaneous decoupling, is given by 



h _ j£ + 2)\ 



dk 



dx 



dTk ji{x) 
dx x 2 



(104) 



where this time Tk is the transfer function of the tensor 
modes. Up to a factor ~ (£/x) 2 this is just what the 
vector modes would contribute if the non-integrated term 
in (94 1 were negligible. In fact, because each power of x 
in the integral over momenta typically yields a factor £, 
this scales with I in the same way as the contribution 
from the integrated term of the vector modes. 

The amplitude of the tensor modes remains constant 
on superhorizon scales and decays also with 1/et inside 
the horizon. Thus, for modes that cross the horizon after 
decoupling but well before the present, the analogue of 
the large-scale approximation (97) is 



dx 



dT k 3e{x) 
dx x 2 



7fe(?7dec) 



M^dec) 



(105) 



Substituting then Eq. (105) into (104) and following the 



same steps as before we obtain 



Comparison with Tensor Modes 



(106) 



It is also illustrative to compare the contribution of 



the vector modes to the angular power spectrum ( 94 ) to 



which, again, holds for £ < 50. Thus, on large scales 
the shape of the spectrum for tensor and vector modes 
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roughly agree if the spectral indices are the same, and 
there is no growing mode in the vector sector. On small 
scales, the situation is different though. Along the same 
lines as before, the contribution of the integrated term 
can be estimated qualitatively by Eq. (105). Carrying 



scales. It follows from Eqs. ( 81 ) , ( 100 ) and ( 101 ) that the 



out the same approximations as for the vector modes, 
the angular power spectrum from the tensor modes then 
becomes 



1(1 + \)C h t 
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(107) 



This result agrees qualitatively well with numerical sim- 
ulations |52l I54| , and lends further support to the ap- 
proximation ( [96] ) for the vector modes on small scales. It 
shows that, at these scales, the angular power spectrum 
of vector modes with spectral index n v is essentially the 
same as the angular power spectrum of tensor modes with 
spectral index 



1 - 2v r . 



(108) 



A particularly relevant example of the equivalence 
arises for a + cu = 0, which leads to v r = 3/2, v m = 5/2. 
In this case, the spectral index of the vector modes is 
n v rj 2, which is just the spectral index of a massless 
field in flat spacetime. Even though the amplitude of this 
spectrum is negligible on large scales at the end of infla- 
tion, because there is a growing mode during radiation 
domination, the amplitude of the spectrum at decoupling 
may be sizable. In any case, during radiation domination 
the spectrum at time is r\ is proportional to (krj) 2 , so all 
modes enter the sound horizon c v kr] = 1 with the same 
amplitude. After horizon crossing the amplitude decays 
as 1/a, as for tensor modes. Hence, up to the frequency 
of the acoustic oscillations, this case is almost equivalent 
to that of tensor modes with a nearly scale-invariant pri- 
mordial spectrum, as stated by (108). The equivalence 
also extends to large angular scales. With v rn = 5/2, 



Eq. ( 101 ) yields a flat plateau in the contribution of vec- 



tor modes to the temperature anisotropics, just as for a 
scale-invariant spectrum of gravitational waves. 



Comparison with Observations 

Present measurements of the CMB temperature 
anisotropies seem to be well-fit by a nearly scale- invariant 
spectrum of scalar perturbations [56] ■ Therefore, if vec- 
tor modes do contribute to the temperature anisotropies 
at observable scales, their contribution must be subdomi- 
nant. This requirement places constraints on the param- 
eters of aether theories, which follow from demanding 



(109) 



Let us obtain a very rough estimate of the contribution 
of vector modes to the temperature anisotropies on large 



bulk part of the contribution stems from the four large 
factors 
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(110) 



Here, we have ignored the (recent) stage of cosmic accel- 
eration, the difference between equality and decoupling, 
and the redshift to the time of equality of matter and 
radiation densities. Using Eqs. (l79| and (81b) we find 



2u T 



1 



-4 + y/l + An, + 3\A + (4/9)«v rj -K r 

_ (HI) 
where k is defined in ( 43 1 , and the indices i and r refer to 



inflation and the radiation era respectively. We have also 
expanded for small k in the last step, and neglected Kj in 
front of n r , because of a relative slow roll suppression fac- 
tor. The sign of k is determined by the sign of (l+a/ci4). 
Hence, if (1 + a/cn) > 0, vector modes are primordially 
suppressed, and the subsequent growth during radiation 
domination cannot compensate for this suppression. On 
the other hand if (l + a/cu) < 0, the growth during radi- 
ation domination may bring the signal well above what is 
observationally allowed. Hence, it seems that the range 
which is the most interesting from the point of view of 
observation is when |1 + a/cu\ -C 1, which corresponds 
to k <C 1. Therefore, it is not excluded that the present 
amplitude of these modes is quite sizable, producing de- 
tectable signals in the CMB, or dipole contributions to 
the gravitational potentials of massive bodies through 
the effect of the velocity field of the vector modes of the 
aether with respect to matter, as discussed at the end of 
the previous section. 

An interesting question is whether, in the range where 
|1 + a/en] <C 1, the contribution of the scalar isocurva- 
ture mode to observables will be larger or smaller than 
that of the vector modes. As we noted around Eq. pHl), 



the relative amplitude of the longitudinal to the trans- 
verse velocity field power spectra is of order ci4/(l + to), 
which is likely to be quite small if the scale of Lorentz 
symmetry breaking is low. Nonetheless, depending on 
parameters, both situations seem possible. Furthermore, 
in a theory such as BPSH, the vector mode is completely 
absent, and we only have the scalar contribution. A full 
analysis of the CMB signatures for the scalar mode is left 
for further research. 



VII. SUMMARY AND CONCLUSIONS 

In this article we have studied cosmological pertur- 
bations in Einstein-aether theory, where the scalar and 
transverse vector sectors of general relativity arc enlarged 
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by an additional dynamical field each. We find that infla- 
tion can induce sizable perturbations in both of these new 
massless fields on observable scales. Our analysis also ap- 
plies to the low energy limit of BPSH gravity, where the 
transverse vector is missing by construction |28) . 

For the purposes of summarizing our results, we shall 
assume that the aether parameters Cj (i = 1, . . . , 4) are 
small. This is natural, since they can be thought of as 
proportional to the square of the ratio of the symmetry 
breaking scale M to the Planck scale c, ~ (M/Mp) 2 <C 1. 

To motivate the choice of the range of parameters we 
shall use below, let us recall that in the Einstein-aether 
theory, the effective gravitational constant on small scales 
Gjv can be different from the effective gravitational con- 
stant which appears in the Friedmann equation G cos . We 
shall call a and C14 the parameters which relate these 
two constants to the bare Newton's constant G. They 
are given in Eqs. ([9| as linear combinations of the stan- 
dard Ci. In terms of a and C14 the effective gravitational 
constants are given by 



G = (l - |) G cos = (l + 



N- 



(112) 



Note that for a + C14 — we have G cos = Gjv- The 
difference G cos — Gat is constrained by nucleosynthesis to 
be less than 10 %, so it seems natural to consider the 
range 



|k| <C 1, where k 



1 



a 



(113) 



This range guarantees the similarity of G cos and Gat, 
but it is typically more restrictive than required by the 
nucleosynthesis bound, since C14 ~ (M/M P ) 2 is naturally 
small. If the parameters are such that we are outside of 
the range ( 113 ), the effects we are investigating would be 



cither too small to be of phenomenological interest, or 
too large to be compatible with observations. 

The main results of the paper are the following. First, 
we find that in the scalar sector, aside from the standard 
adiabatic mode £ (which corresponds to the curvature 
of surfaces of constant matter density), there is an ad- 
ditional isocurvature mode which can be important for 
phenomenology. Geometrically, the isocurvature mode 
can be described as the differential e-folding number SN 
which separates the surfaces of constant matter density 
from the surfaces orthogonal to the aether. This plays the 
role of a velocity potential v for the aether with respect 
to matter. At the time of horizon exit during inflation, 
the amplitudes of SN and v are comparable to that of 
the standard adiabatic mode £: 
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Mp 



-1/2 



(horizon exit). 



(114) 



Here H is the Hubble rate and e <C 1 is the slow roll 
parameter during inflation, which is independent of the 
aether parameters. 

After horizon crossing, the curvature perturbation £ 
stays constant, while the behaviour of SN depends on 



the parameter k defined above. For k < 0, the isocur- 
vature perturbation slowly decays on large scales, while 
for k > it grows. On the other hand, the velocity 
perturbation is given by v ~ (k/a)SN, where k is the co- 
moving wave number and d is the derivative of the scale 
factor with respect to proper time. Hence, during infla- 
tion, when d grows, the long wavelength velocity field 
decays, roughly in proportion to the inverse of the scale 
factor. After inflaton, the universe decelerates and the 
velocity field grows again. At the time of horizon reentry, 
on cosmologically relevant scales, we have 

v ~ SN ~ e^C ~ e Nk 1CT 5 < 1. (horizon reentry) 

(115) 

where N ~ 60 is the number of e-foldings of inflation 
since the time when the cosmological scale first crossed 
the horizon. The last inequality indicates the limit of 
validity of the linear approximation. Note that for k = 0, 
the isocurvature perturbation and the velocity field of the 
aether are comparable to £ ~ 10~ 5 at horizon reentry. 
However, with k < 10/A, we can have SN < 1. If k is 
large enough to saturate the inequality, this still allows 
for mildly relativistic speeds for the aether field v ~ 1 
within the observable universe. 

Similar results hold for the vector sector. Denoting 
by V the transverse component of the aether velocity 
field with respect to matter, we find that on superhorizon 
scales 



V 



VC14 



1/2 



Hence, if C14 < e (which seems quite natural if the scale of 
Lorentz symmetry breaking is low), the vector contribu- 
tion to the velocity field will be dominant with respect to 
that of the longitudinal component. On the other hand, 
in a theory such as BPSH, the transverse component V 
is missing, and the scalar part v is the dominant one. 

We also find that the longitudinal gauge gravitational 
potentials <f> and ip can be different even for the adia- 
batic mode. On superhorizon scales, we find that this 
effect (which can be attributed to anisotropic stress of 
the aether energy momentum tensor) is of order: 

(</> ~ VOadiab ~ ^adiab C 13 ~ £ C 13 ~ 10" 5 C 13 . (116) 



where C13 ~ (M/Mp) 2 is another combination of the 
aether parameters Cj, given in Eqs. Physically, this 
parameter can be expressed in terms of the propagation 
speed of tensor modes C13 = c^ 2 — 1. The isocurvature 
mode contributes maximally to the anisotropic stress, 
but the potential due to the isocurvature mode is sup- 
pressed by C13: 



C13 SN. 



(117) 



Since SN can be larger than (, the anisotropic stress can 
be dominated by the isocurvature mode. The anisotropic 
stress on observable scales is suppressed from its value at 
horizon crossing, due to the dynamics of the aether on 
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subhorizon scales. For it = 0, the effect scales like k~ 2 
for modes that crossed the horizon during the matter era. 
For modes that crossed the horizon during the radiation 
era, the behaviour changes to A: -1 . Current constraints 
on <f) — ip on cosmological scales are not very restrictive, 
and | C13 1 < 1 seems to be allowed by observations. 

The aether manifests itself in PPN parameters through 
frame dependent effects, which cause anisotropies in the 
gravitational field of bodies which move with respect to 
the aether. In this way, the velocity field generated dur- 
ing inflation might be detectable. It should be noted, 
however, that it seems difficult with present technology 
to observe the statistical properties of the random field 
from this particular type of observations. Even if the ve- 
locity field were relativistic on cosmological scales v ~ 1, 
it falls with scale as k~ 2 . In particular, the component 
which varies on scales of the order of 100 Mpc would then 
be below the virial velocity v v - lr ~ 10~ 3 of objects bound 
in galaxies, and it seems unlikely that we can directly 
sample frame dependent effects in objects which are lo- 
cated at distances larger than that. On the other hand, 
at the relatively small distances where the observation 
of frame dependent effects is accessible, we may still de- 
tect a large but fairly homogeneous velocity field, even 
one that is much larger than the virial velocity of bound 
objects. 

Finally, we have computed the contribution of trans- 
verse vector fields V to the angular power spectrum of 
CMB anisotropies. We find that for k = 0, the spectrum 
of multipole coefficients Cj has the same shape as that 
of tensor modes. The amplitude, on the other hand, is 
related to the spectrum C\ for tensor modes and C\ for 
the adiabatic scalar mode as 



C) 



c 13 Q 2Nk 
C14 



C) 



-13 n 2Nkr<(, 



£'14 



(118) 



This means that the vector modes in Einstein-aether the- 
ory can easily dominate the signal from tensor modes. 
The analysis of polarization induced by the vector modes 
is therefore of phenomenological interest, and is left for 
further research. Moreover, we know that the CMB is 
well-fit with a primordial spectrum of scalar adiabatic 
perturbations. This imposes additional phenomenolog- 
ical restriction amongst the parameters C13 and R of 
Einstein-aether theories, of the form 



1 

2N 



In 



Cl4 



"13 



(119) 



So far, we have not included the constraints which fol- 
low from the frame-dependent effects on the PPN param- 
eters. These are summarized in Appendix |Aj and take 
the form 



wai<10" 7 , oj 2 a 2 < 10" 



13 



(120) 



Here, to = max{V, v, v v - lT }, is the velocity of the aether 
with respect to the object whose gravitational field is 



the typical virial velocity for bound objects with respect 
to the CMB frame. The post-Newtonian parameters a\ 
and a 2 are combinations of the four aether parameters 
(a, C14, c+, c_). Here, following [57] . we have introduced 
c + = C13 = C1+C3 and c_ = c\ — C3. Phenomenologically, 
it is possible to set a\ — a 2 = 0, which determines a and 
C14 as functions of the other two parameters in the model, 



a = — C14 = —2 



{C++C-Y 



(121) 



The parameters C+ and c_ remain rather unconstrained 
by observations. Stability requirements and superlumi- 
nality (or Cherenkov) constraints are satisfied provided 
that — 1 < c + < 0, c + /3(l + c+) < c_ < 0. Constraints 
from radiation damping in binary systems determine fur- 
ther constraints on the (c_|_,c_) plane, but a sizable co- 
efficient 



C13 



< 1 



(122) 



still seems to be allowed by all observations [S7]. This 
is important, since the gravitational effects of the aether 
are suppressed by this coefficient. For instance the con- 
tribution of vectors to the CMB anisotropics is of order 



C 



r 2 V 2 
c 13 v J 



(123) 



where V < 1 is the aether velocity field. Hence, the 
observability of the effect depends crucially on C13 being 
sufficiently large. 

This brings us to the question of fine tuning amongst 
the parameters of the model. In a low energy theory, one 
might have expected all dimensionless parameters to be 
of the same order, 

a ~ (M/Mp) 2 . 

Observability of Cj requires an inequality of the form 
Ci3^ <L 10~ 6 , which would be natural provided that 



(M/M 



Pjobs 



> 10~ 6 W 



(124) 



On the other hand, in Eq. (121) we have adjusted the 
para mete rs so that u\ = a 2 = 0, but the actual restric- 
tion (120 1 is of the form a 2 ^ 10 -13 w -2 . Hence, a 2 must 



be well below the natural scale (124) by a considerable 
suppression factor 



a 2 < lO^w-^V/a;) {M/Mp) 2 ohs , 



(125) 



being tested at post-Newtonian order, and v v 



10- 



with 10 -3 < lu < 1. In the classical theory, the parameter 
a 2 can always be chosen by hand to have any particu- 
lar value. However, in an effective field theory (EFT) a 
parameter is considered to be finely tuned or technically 
unnatural if quantum corrections to it are larger than 
the desired renormalized value of the parameter. The 
question, therefore is whether the very small values of 
a 2 -C (M/Mp) 2 are stable or not under quantum correc- 
tions. Withers [SS] has recently analyzed the Einstein- 
Aether theory as an EFT. with the conclusion that the 
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parameters Ci receive only negligible logarithmic correc- 
tions. A similar result may hold in BPSH theory [27] . 
This subject is left for further study. 

To conclude, the results presented here show that the 
preferred frame singled out by the aether field A 11 , or 
by the preferred foliation of the BPSH theory, may have 
picked up a large random velocity field seeded by quan- 
tum fluctuations during inflation. Depending on the pa- 
rameters, this may even be mildly relativistic on cosmo- 
logical scales. The effects of this velocity field may be 
detectable in observations of frame dependent PPN ef- 
fects, or in specific features in the CMB spectrum such 
as a sizable contribution from vector modes. These issues 
deserve further investigation. 
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Appendix A: Other Constraints 

In this article we have derived the constraints on the 
parameters of Einstein-aether theories that follow from 
classical and quantum stability, and from phenomeno- 
logical considerations related to the primordial perturba- 
tions. In addition to these constraints, there are further 
conditions that the Cj's have to satisfy, which arise from 
the Post-Newtonian limit of the theory, Big-Bang nu- 
cleosynthesis, and from the arrival of high-energy cosmic 
rays to earth. An extensive summary of these constraints 
can be found in [57] , 



Post-Newtonian Limits 

In any metric theory, the gravitational field created by 
non-relativistic bodies can be characterized beyond the 
Newtonian limit by a set of post-Newtonian PPN pa- 
rameters, whose values are tightly constrained by solar 
system tests of gravity [5] . The parameters j3 and Cxa in 
aether theories agree with those of general relativity, and 
also agree with the measured ones [I3J . But because the 
aether defines a preferred frame, it also introduces ad- 
ditional departures from general relativity, which mani- 
fest themselves as gravitational potentials that depend on 
the velocity of the interacting bodies with respect to the 
aether. These preferred-frame effects are encoded in the 



PPN parameters a.\ and a,i- One of the most stringent 
limits on the value of ai comes from measurements of 
the eccentricity of the binary pulsar J23I7+I439 (which 
would change if a\ were non-zero) [58] . while one of the 
most stringent limits on a-i stems from the alignment of 
the sun spin with the solar system angular momentum (a 
non-zero a-i would lead to a misalignment) |59) . These 
limits lead to the conditions 



a 2 = 



(ClC 4 



4) 



2ci 
ai _ 
2 



< 1.7 x 10" 



(Ala) 



(2C13 ~ C14)(Q + Cl4) <1.2xl0- 7 . (Alb) 
p{2 - ci 4 ) 



It is important to stress that both limits assume that the 
velocity of the sun with respect to the aether uj is the 
velocity with respect to the frame in which CMB dipole 
vanishes, ui ~ 10~ 3 . Roughly speaking, the limit on ot\ 
is actually a limit on a\UJ, while the limit on o?2 actually 
constraints o^w 2 . Thus, if lu were larger that assumed, 
as the results of our work seem to allow, the limits on 
«i and a.2 would be actually tighter. In other words, the 
constraints on the PPN parameters a\ and ai actually 



Q-x < KT 4 



io- 



Q' 2 



< 10" 



id 

10- 

ui 



(A2a) 
(A2b) 



where lo is again the velocity of the sun with respect to 
the preferred frame. The constraints (All are typically 
satisfied if the norm of the A, defined in Section [n] is of 
order M ~ 10~ 4 M P . 



Big-Bang Nucleosynthesis 

The agreement between the predicted light element 
abundances, and those actually observed (or indirectly 
measured using CMB observations [56]) constrains the 
value of the Hubble constant at the time the light el- 
ements formed, at temperatures of about T w 10 9 K. 
Because the expansion rate depends on the value of 
the renormalized Newton constant G cos = 2G/(2 — a) 
through Eq. (12a), and because the latter is related 
to the "Newtonian" gravitational constant by Eq. ( |14[ ), 
given the measured value of Gn on small scales and the 



number of relativistic species during nucleosynthesis, one 
can determine how light element abundances depend on 
the parameters a and ci 4 . Agreement of such a predic- 
tion with observations then implies 



Cl 4 



- a 



< 10%. 



(A3) 
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Cherenkov Radiation 

If any of the propagation speeds of tensor, vector or 
scalar modes we have discussed were sufficiently smaller 
than the speed of light, highly relativistic particles trav- 
eling close to the speed of light would loose energy into 
these modes by a process analogous to Cherenkov radia- 
tion (this is kincmatically possible only if the dispersion 
relation of the emitted quanta is subluminal.) Although 
the emission amplitude is inversely proportional to the 
Planck mass, the fact that we detect these cosmic rays, 
and that they must originate at astrophysical distances, 
allows one to place quite stringent limits on the parame- 
ters of the aether l60l. 



C13 < 1 x 10 
< 1.4 x 10 



lo 



c 13( c 13 + 2c 4 ) f i m ;.| 



(c3 ,-? )2 < 1 x 10- 30 



|C14| 
C4 - C2 - C 3 
Cl 



< 3 x 10 
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(A4a) 
(A4b) 

(A4c) 
(A4d) 



It is important to realize though that these constraints 
only apply if the different aether modes pro paga te sublu- 
minally. Under this assumption conditions ( A4 1 can also 
be taken to imply the bound M < 10~ 7 Mp on the norm 



of the aether field A M defined in Section |II| 



Propagation Speed 

Some authors impose further conditions on the param- 
eters of the aether, namely, that the propagation speed 
of the perturbations be subluminal. The origin of this 
requirement goes back to the violations of causality that 
appear in Lorentz-invariant theories with superluminal 
signals. As far as we know, there is no link however be- 
tween superluminality and violations of causality in back- 
grounds like the ones we are considering. The cosmic 
aether breaks Lorentz invariance and defines a preferred 
reference frame. Signals always travel forward in time in 
this frame, so no closed timelike curves can arise. Even 
the construction of |61j , in which due to the nature of the 
background closed timelike curves may appear seems dif- 
ficult to realize here, because the aether satisfies a fixed- 
norm constraint. Hence, we shall not require subluminal 
propagation, though because this is a somewhat contro- 
versial issue, we collect the appropriate conditions here 
for completeness. They easily follow from Eqs. (27), ( p7| ) 
and (74). In the limit of small coefficients, Ci <§; 1 they 
read 



ci - c 4 - 



C13 > 



i 3 > o 

/3-ci4 > 0. 



(A5a) 
(A5b) 
(A5c) 



Note that because we have taken metric perturbations 
into account, these conditions differ from those derived 
in the limit of a decoupled aether [33]. For alternative 
views on superluminal propagation we refer the reader to 
the references 



Appendix B: Scalar equations in the longitudinal 
gauge 

In this appendix we discuss the scalar sector in the 
longitudinal gauge. This is useful and complementary to 
the next appendix, which deals with the gauge-invariant 
formulation. 

In the longitudinal gauge, there are five fields in the 
scalar sector: 4>, ip, C, SX and the inflaton perturbation 
Sip. Therefore, we need five independent equations to 
uniquely determine their values. Contracting Eq. ^ 
with A@ and using the constraint A 13 Ap = — 1 yields an 
equation that expresses the Lagrange multiplier in terms 
of the aether and the metric, 



SX 



1 



6(a - c 2 )U 2 



6c 2 



3c 2 H(f>'- 



+ c 3 k 2 cf> - 3(2^ - c 2 )Hip' + 3c 2 V>"- 
- OS + c x )Uk 2 C +((3- ci)fc 2 C 



, (Bl) 



which shows explicitly how the Lagrange multiplier can 
be expressed in terms of the remaining fields. The time 
component of the linearized aether field Eq. ^ is iden- 
tically satisfied. The linearized spatial components give 



C" + 2HC' - 

-k 2 C 



a 



C14 



2W 



a 
a 



a 
a 



C+ 



Cu 



H(j> + <j>' + — ip' = 0, (B2) 

Cu/ C14 



which combined with the ° i Einstein equation results in 







2 - a 



k 2 C = 4ttG ip'S(p. (B3) 

2 — a 



Eq. (B3) expresses Sip in terms of the remaining scalars, 



and allows us to eliminate Sip from our system of equa- 
tions. On large scales, this equation has the same form it 
would have in the absence of the aether, with the differ- 
ence that the effective Newton's constant has the renor- 
malized value implied by Eqs. (12). The part of the % 



Einstein equations which is not proportional to S % , is 



<t> = ip + c 13 (C' + 2HC), 



(B4) 



which immediately reveals that the Einstein-aether is a 
source of anisotropic stress in the scalar sector. This 
equation allows us to express 4> i n terms of tp and C, and 
thus eliminate yet another variable from the equations. 
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Note that scalar fields and perfect fluids cannot source 
anisotropic stress, which is why a value of ip — <p different 
from zero is sometimes attributed to modified gravity. 
Finally, the sum of the % and the l j Einstein equations 
proportional to S l j is 



2- a 



2- a 



k 2 tp- 



2 - a 



-k z C 

8ttG 

2- a 



2- a 



Uk z C 



m(l-w)<p'5ip, (B5) 



where we have used Eq. (Bl| to eliminate <5A, and that 



during power-law expansion the equation of state param- 
eter w is constant. 

Eqs. dBlL (lB2j), (|B3l, (TbII) and (|B5]) form a set of 



five differential equations for the five unknowns. We can 
use the constraints ( B3 ) and ( B4 ) to eliminate <p> and 5ip 



from Eqs. (B2) and (B5), arriving at 



C" 



and 



2+ C ^ + a } )uC>- 

Cl4(l + C13) / 

/ ci4 - a + 2c 13 ci4 a" 
\ c 14 (l + c 13 ) a 



2a, 
cu 



n 2 c+ 



P ,2 r < , c 14 + a 



Cl4(l + C13) 



k z C + 



c i4 (l + C13) 



(B6) 



1. Short-wavelength Solutions 

In the short-wavelength regime, k\ri\ 3> 1, the solutions 
of the equations of motion (B6) and (B7) behave approx- 



imately like in flat space. The notion of an approximate 
solution can be formalized by introducing krj as an ex- 
pansion parameter. In the limit k\rj\ 3> 1 the solution of 
Eqs. (B6) and (B7) can be cast in the form 



ip = ij(kri) exp(-ic s krj), C = \c{kr]) ip, (B8) 

where, ip and C are functions whose power series expan- 
sion starts at a finite positive power of krj, and c s is a 
"sound speed" to be determined. Substituting the ansatz 
(|B8|) into Eqs. (|B6|) and (|B7|), and keeping the leading 



powers of kr\ yields a set of algebraic equations with two 
positive frequency and two negative frequency solutions, 
for a total of four solutions, as expected. At leading or- 
der, ip remains unconstrained and can be taken to be 
constant. The positive frequency solutions are given by 

2 



(C s )a : 



C a = i 



a 



P 



P - Cu(l + C13) ' 



(B9a) 
(B9b) 



where c a is the sound speed of Eq. ( 37 1 



These two modes correspond to the two independent 



short wavelength solutions ( 39a I and ( 39b I that we found 



in Subsection |IV A| To see that this is the case, we may 
use the expression of 5N and £ Q in the longitudinal gauge 



ip" + 3(1 + w)Hip' + 2 2H 2 )ip + 3(1 + w)H 2 ip+ 



c 13 HC" + 2c 13 H 2 C" + 3(1 + w)c 13 U 2 C'+ 



C13 ( 6^- - ion 2 ) nc + 6(1 + w) Cl3 n 3 c+ 

2 + Cl4 ;„2„,. , c 14(l + Ci 3 )~/3 fc2c ,/_ 



2 - a 



■kfip- 



2 - a 



c 14 (l + 2c 13 ) + 4/3 - a - 3(1 + w)/3^ l2 



2- a 



Hk z C = 0. 



(B7) 



Because this is a system of two second order linear dif- 
ferential equations, we need to specify four independent 
initial conditions, so there must exist four linearly inde- 
pendent solutions. This is also what we expect by simply 
counting matter fields. In the limit of weak gravitational 
couplings, we may neglect metric perturbations, so we 
just have one degree of freedom in the inflaton perturba- 
tions and one degree of freedom in the aether field pertur- 
bations, for a total of four initial conditions to determine 
uniquely the evolution of the system. As we deviate from 
the limit of weak coupling, neglecting metric perturba- 
tions ceases to be a good approximation, but the number 
of degrees of freedom in the theory remains unchanged. 



SN 



HS(p 



(a = $-UC. 



nc, 



(BlOa) 
(BlOb) 



In the first equation, dip should be expressed in terms of 
ip and C through the relation 

M 2 

Sip = ^ [(2 - a) (iP' + c 13 UC + 2c 13 n 2 C) - /3k 2 

(Bll) 



which follows from ( B3 1 and ( B4 ) . By comparison with 



( 39a I and ( 39b I we also obtain the overall normalization 



factor ip. For the first mode, we have 

Ipa 



Zl/2 —ic a krt 1 
a C 1 ~ 

a y/2c a k k 



(B12) 



where Z a is given in Eq. ( 36 1 , and C a in Eq. ( B9a 
the second mode, we have 



For 



^■1/2 —ikr] ^ ^ 



2k 



where C v is given in Eq. ( B9b I and 



Z 1/2 = i- 
* ~ kMl 



Cui/3 + a{l + c 13 )) 
P - Cl 4 (l + c 13 ) 



(B13) 



(B14) 
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Substitution of (B12) into (BIO I reproduces Eq. (39b I, 



while substitution of (|B13 1 into (BIO), together with the 
background Eqs. (fl2~E reproduces Eq. (39a). The vac- 



and the difference of the two metric potentials (which is 
proportional to the anisotropic stress) are given by 



uum is thus characterized by the two independent solu- 
tions of Eqs. (B6) and (B7) that approach (B12) and 



(B13) in the limit k\rj\ — > oo. 



. (5 + 3ttf)(l + ci 8 ) . 
Ci = wrr-, r-; tz , » J i, 



C2 



3(1 + to) 
0, 



ci 3 (5 + 3w) 



fl>i = - 



Cl3 



<1 



I + C13 

(B19a) 
fa = 0. (B19b) 



2. Long- wavelength Solutions 

In the limit of long wavelengths, k\rj\ -C 1, we ma y 
neglect terms proportional to k 2 in Eqs. (B6) and (B7). 
In this limit, the power-law ansatz 



*l> = (-»?)*, C = C-(-ri)-1) 



(B15) 



reduces the two coupled differential equations ( B6 ) and 



(B7) to an algebraic system for the two constants t and 



It can be readily checked that for these modes SN = 0, so 
that matter is at rest in the aether frame. Though these 
adiabatic modes have the properties described in [39] . 
they do not share the properties postulated in |T0TI72*] . In 
particular, for the first adiabatic mode, the anisotropic 
stress is non-zero. The form of the two adiabatic modes 
for an arbitrary expansion history and matter content is 
derived in Appendix |P| 



b. Isocurvature Modes (C, = 0, SN 7^ 0) 



2 . cu(l + ci 3 )(5 + 3w) + 2ci 3 (ci4 + a) 

t ~r ~ m ~ — ; 1~ 



+ 



c u (l + c 13 )(l + 3w) 
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Cl 4 (l +C13) 
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2 

l + 3w 



(B16a) 



The two remaining solutions of Eqs. (B16) require 
C = (1 + 3u>)/2ci 3 , which gives 

From (B16a), the exponents are given by 



(B20) 



t(5 + t + 3w + 3wt)(l + 3w - 2Cci 3 ) = 0. (B16b) 




4k, (B21) 



2t+ = - 



where k is given by Eq. ( 43 ) . It is straightforward to 



Because Eqs. (B16) arc linear in C, they may be reduced c h ec k that for these modes we have 



to a single quartic equation for t, with four different so- 
lutions, as it should be. 



Adiabatic Modes (5N = 0) 



C± = 0, 

; - ip± = [ 1 \ 3w ) t {T ) i>± 



(B22a) 
(B22b) 



SN+ = - 



Two solutions of the coupled equations (B16) follow 



directly from Eq. (B16b), 
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U =0, 
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1 + 3w 



to = ~- 



h + 3w 
l + 3w' 



c 2 = - 



3(l + w) + ci 3 (5 + 3u;) : 
1 + 3w 



(B17a) 
(B17b) 



These are two isocurvature modes, in the sense that the 
curvature perturbation on comoving slices £ vanishes for 
any value of w, 

From Eq. (B21) it is straightforward to check that, 



The corresponding perturbations are the two "adiabatic" 
modes that always exist at long wavelengths, regardless 
of the matter content of the universe [39]. Along these 
two modes, the (spatial) curvature perturbation on co- 
moving slices, 9 



C = ^ + 



3H{l + w)' 



(B18) 



for any value of w, one of the two modes is a decaying 
one. Whether the second mode is growing or decaying 
depends on the sign of k, which is in turn determined by 
the sign of 1 4- (a/cu). For k > the second solution 
is also a decaying one, but for k < there is a growing 
mode. In the special case k = 0, there is a constant 
non-decaying long wavelength solution. 

The existence of a growing non-adiabatic isocurvature 
mode in Einstein-aether theories for (a/cu) < —1 can 
have important phenomenological consequences, as we 
discuss in the main text. 



Recall from equation | [32| that we mean comoving with respect 
to all forms of matter, excluding the aether. The °j Einstein 
equation | |B3| however reveals that the contribution of the aether 
to the total velocity perturbation is negligible on large scales. 
Hence, on large scales, hypersurfaces comoving with matter and 
comoving with matter plus aether are actually the same. 



Appendix C: Canonical reduction of the scalar 
sector. 

The normalization of the spectrum of scalar pertur- 
bations follows from the normalization of the action for 
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the corresponding physical degrees of freedom. Here, we 
find the reduced set of gauge-invariant dynamical vari- 
ables, and express the second order Lagrangian in terms 
of these. This Lagrangian can also be used, of course, to 
rederive the scalar equations of motion (B6| and (B7|. 



The starting point is the Lagrangian for scalar per- 
turbations in an arbitrary gauge, which is obtained by 
substituting the metric plj ) into the action ([3]), and ex- 
panding to second order in the scalar perturbations. Us- 
ing the constraint ^ which is obtained from variation 
with respect to 5X, we arrive at 



£(2 ) = ^P fl2 



Ml 
2 



2k 2 ^ 2 - 3(2 - a)?// 2 - 4fcV</> + 4k 2 iP'B- 

1 



- (2 - a)k 2 ^'E' + 2ak 2 ip'C + /3k 4 {C + -E') 2 

- c 1A k 2 {cf) + C + B'f - 6(2 - a)H(M/- 

- 2(ci4 - 2)Hk 2 cj)B - (2 - a)Uk 2 4>E'+ 

+ 2(a - c 14 )Hk 2 (/)C - (2 - a)(2H 2 + H')(f> 2 + 
+ {a(H 2 - W) + c u (H 2 + n'))k 2 {C + B) 2 + 



+ M p 2 ( <V 2 - k 2 Sip 2 + 2ip'6<p(3<ip' - k 2 B + 



k 2 E 



2ip'5ip'<j> ~ o J V w V - 2a 2 V ip 6<p<j>) 



(CI) 



Not all variables in this Lagrangian are dynamical. Some 
linear combinations are gauge modes, while others are 
constrained. We would like to find a Lagrangian that 
contains dynamical gauge-invariant variables only. 

The identification of constraints and the reduction of 
phase space is best performed in the canonical formalism, 
where the equations of motion are at most of first order in 
time. Constraints are equations of motion without any 
time derivatives, and can be substituted back into the 
first order Lagrangian. Here, we closely follow Fadeev 
and Jackiw's method for dealing with constrained sys- 
tems [64]. For a discussion of cosmological perturbation 
theory in this framework, see |65) . 

We begin by introducing new variables U and W 
through 



2W = B + C, 
2U = B-C. 



(C2a) 
(C2b) 



The conjugate momenta of the system are given by 



(2) 



n. 



c_ 



Mia 2 



+ 3M p 2 ip'S<p 



ak 2 (W -U) + 2k 2 (U + W)- 

k 2 E' -3(2-a)(V/ + H<p)+ 
(C3a) 



( 2 " a ) ,,.-> w 



n, 



BE' 



CP 
dSip' 



1 



-Mla 2 k 2 



Pk 2 (W -U + -E')- 



- (2 - a)(ip' + H<j>) + Mp 2 Lp'5ip 



a 2 {8<p' -fa'), 



(C3b) 
(C3c) 



C 



(2) 



dW 



- 2c u M 2 P a 2 k 2 ((t> + 2W'), (C3d) 
and we can write the first order Lagrangian 

cm = u e e' + u w w' + n AV <y + n. v'- 

n 2 



3U% 



n 2 



+ 



+ 



M 2 a 2 fc 4 (l + ci 3 ) 8ci 4 fc 2 M 2 a 2 2a 2 

2k 2 f3WU i> 

4M 2 L a 2 (2 - a)(l + c 13 ) ~ (2 - a)(l + c 13 ) 

2(1 - c 13 )wn E 

(l + cis) 



A/ 2 a 2 fc 2 (l + c 13 ) 



2UU E 



M 2 (2-a) 



+ Mpa 2 k 2 i> 2 + 2Mpa 2 k 2 W 2 x 



2k 2 (3 



\ v (2-a)(l + c 13 ) 

\ aHlp2 ( fc2 + ° 2 ^- M 2 (2-a) 
2aa 2 k 2 ip , 5ipW 



(c M + a)H 2 + (c 14 - a)H' 
3</ 2 



2c u Mpa 2 k 2 HW 



(2 -a) 



1 



- 2Mpa 2 k 2 tp + -U w + HH^, - <p'ILs v 
-a 2 5Lp{a 2 V v + mip')}. 



(C4) 



Variation with respect to the independent variables ip, 
n.0, E, He, 5<p, lisp, W, Hw, U and <f> leads to the same 
equations of motion as those derived from the variation 



of (CI) with respect to ip, E, dip, B, C and 



Note that the time derivatives of U and <\> do not ap- 



pear in Eq. ( C4 ) , so variation with respect to these vari- 



ables leads to the two constraints 

n s = o, 



(C5a) 



n. 



Sip 



Ac 1A M 2 a 2 k 2 HW - 2a A V v 5ip - AM 2 p a 2 k 2 i) 



2ip' 

U w + 2-HH0 - 6a 2 H(p'8ip 



(C5b) 
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Substitution of these constraints also causes E, <p and U 
to drop from the Lagrangian, which therefore depends 
only on the five independent canonical variables ip, H^, 
Sip, W, Hw- Five is one too many, since we expect two 
canonical pairs only. Indeed, one of the variables is re- 
dundant, and it corresponds to the residual gauge in- 
variance of the Lagrangian. Let us introduce the gauge- 
invariant combinations 



sn = rm = 2hw + —Sip. 



(C6a) 
(C6b) 



Geometrically, these can be interpreted as follows. The 
variable (, is the curvature perturbation on surfaces of 
constant inflaton field ip. The variable SN is the same 



as the one introduced in ( 33 ) , and can be interpreted as 



the differential c-folding number between hypersurfaces 
of constant inflaton field and surfaces orthogonal to the 
aether field. 

The momenta conjugate to the gauge-invariant vari- 
ables C and n are given by 



„ „ 2M 2 a 2 k 2 . 

n c = n V; + — — Sep, 



(C7a) 



2c li M 2 a 2 k 2 'H . rr .„ , 

n n = U w + p , Sip. C7b 

ip* 

In terms of the new variables, the field perturbation Sip 
disappears from the Lagrangian ( C4 ) and we have 



CP GI = M 2 a 2 k 2 



c 2 + 



2 I 2 !*!f + (c 14 + a)H 2 + (cl4 - a)W ) Q 2 
i (2-a)(l + ci 3 ) 



2M 2 k 2 
ip' 2 



P 



i 



8c 14 M 2 a 2 k 2 8a 2 <p 

n 2 \ h 



4M 2 a 2 (2-a)(l + c 13 ) 2aV 2 J ^ 2a 2 ^ 2UcUn 



c 14 M P Hk 2 



IP 



(2^n c + n n )fi 



2f3k 2 



2M 2 Uk 2 



ip'- 



M 2 k 2 



(2-a)(l + d 3 ) 

^cnn + n f c' + nnO'. 



(C8) 



Expression ( C8 ) gives the first order Lagrangian we have 



been looking for, since it is a function of two canonical 
pairs, corresponding to two field degrees of freedom. 

To see this more explicitly, we may vary with respect 
to IIq and IT, and plug the resulting equations back 
into Eq. (C8| to obtain the second order Lagrangian. 



For reference we just reproduce the leading terms in the 
limit k\rj\ 3> 1 (the full expression is cumbersome and not 
very illuminating). In terms of £ and SN, this is given 



by 

C (2)GI 



M 2 



-4(2 + c 14 )fc ; 

C14 



asNy 



fc 2 (2(4 + cua) + c 14 (a - 2)(1 + 3w)) 



2c 



14 



2(2 + c 14 )fc 2 / . 2 , 4(2-a)(l + c 13 ) 



'■'14 



P 



(SN)'C 



2(2-a)(l + c 13 ) /2 



(2 - a)(4(l + c 13 ) - 3/3(1 + «;)) ,_ ArV2 



2/3 



-(57V)' 



(C9) 



where the ellipsis denote terms which are subleading in 
the momentum expansion. 

Variation of ( |C9[ ) with respect to C and <5iV (including 
the terms that we do not explicitly write down) yields two 
second order differential equations for ( and SN. These 
equations of motion are valid in any gauge. To find their 
form in the longitudinal gauge, we may use Eq. (BID 



to express the inflaton perturbations in terms of met- 
ric and aether perturbations. Substituting in (C6a| and 



(C6b), we can cast the equations of motion for ( and fl 



as two third order differential equations for the longitu- 
dinal gauge variables tp and C . The latter happen to be 
precisely linear co mbin ations of Eqs. (B6), (B7) and the 
time derivative of (B7|. 10 



In Eq. (C9), the curvature perturbation £ on surfaces 



of constant inflaton field is coupled to the variable SN. 
However, if we replace £ by the curvature perturbation 
(a = ( — SN on hypersurfaces orthogonal to the aether, 
this leads to a Lagrangian for two decoupled variables, 
Cq and SN: 



1 



1 



ZZn ££a 

(CIO) 

"where the ellipsis denote terms which are subleading in 
the momentum expansion, and we have introduced 



_! _ 3(l+iu)(2-a) 



Z N - 



Mia 



pU , 



(Gil) 



and 



_! _ 2(l + ci 3 )(a-2) 



z- L = 



P 



M P a 2 ,ci = - 



(2 + c 14 )/3 



c 14 (a-2)(l + ci 3 )' 
(C12) 

This form of the Lagrangian will be used in order to 
normalize the positive frequency modes associated with 
the initial vacuum fluctuations. 



10 Note in particular that Eqs. | |B6[ | and | |B7[ ) cannot follow from 
a variational principle from a reduced Lagrangian depending 
quadratically on C and if). If a = — C14, the evolution of C 
decouples from that of ip, while the evolution of the latter does 
depend on the evolution of the former. 
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Appendix D: Long Wavelength Adiabatic and 
Isocurvature Modes 



The properties of the two long wavelength adiabatic 
and isocurvature modes for arbitrary expansion history 
and fairly general matter content can be also obtained 
by following a procedure outlined by Weinberg in [39 . 



1. Adiabatic Modes 



Consider the gauge transformations generated by 



r\ — > r\ + e{rf) and x l 



■ uj x 



(Dl) 



where u is a constant. Using the transformation prop- 
erties of the metric one finds that these transformations 
preserve the structure of longitudinal gauge. In partic- 
ular, they induce the following transformations on the 
metric and aether perturbations, 

0-> -He, iJj ^ ip + u + He, C^C + e. (D2) 

Because the equations of motion are invariant under 
gauge transformations, the difference of two sets of per- 
turbations that differ by a gauge transformation is a so- 
lution of the linearized equations, 



He, i/} = u + He, C 



(D3) 



The corresponding values of the remaining perturbation 
variables can be also determined by their transformation 
properties under (Dl). For instance, for any scalar per- 
turbation Sip or any velocity perturbation Sui = didu the 
solutions have 



5ip = —eip / , 6u = ae. 



(D4) 



Of course, these space-independent solutions are just 
gauge modes, physically equivalent to no perturbation at 
all. But they can be extended to actual space-dependent 
perturbations if the linearized n i and \. Einstein equa- 
tions are satisfied for these putative solutions. The 



o . 



equation is automatically satisfied for the ansatz ( D3 ) 
On the other hand 



and (D4) 



aether the 
straint 



in t he presence of the 
j Einstein equation ( B4 1 imposes the con- 



2He 



1 



l + cia 



0. 



(D5) 



where we have assumed that the remaining matter does 
not contribute to the scalar anisotropic stress. The gen- 
eral solution of Eq. (D5) is the superposition of two 
solutions, with 



ei 



<~2 



1 UJ 

a 2 1 + ci 3 



C 



dfja 2 {fi), oji—uj, (D6a) 
oj 2 = 0, (D6b) 



where Co is an integration constant. The first solution 
yields the non-decaying mode, which in the "gravity" sec- 
tor reads 



1 + C13 



C X = -- 



1 - 



H 1 

a 2 1 + ci 3 



dfj a 2 (fj) 
dfja 2 (fi) 



1 



dfja 2 (rj). 



(D7a) 
(D7b) 
(D7c) 



1 + ci 3 a 2 

This reduces to the adiabatic mode ( |B17a I for a constant 
equation of state. For this mode the curvature pertur- 
bation is constant, £ = u>, and the anisotropic stress is 
non-zero (if C13 7^ 0). The second solution in Eq. (D6a) 
corresponds to a decaying mode, which, for a constant 
equation of state, agrees with the adiabatic mode in Eq. 
pl7bl 



= C 



H 



4'2 = C 



H 



C = 



Co 

n 2 ■ 



(D8) 



For this second adiabatic mode, the curvature perturba- 
tion vanishes, £ = 0, and so does the anisotropic stress. 



2. Isocurvature Modes 

An extension of the previous method also unveils the 
two isocurvature modes, under the assumption that the 
aether does not couple to matter. Consider the ansatz 



c 13 {C'+HC), yj = -c 13 HC, 



(D9) 



which arises from the gauge transformation (Dl) with 
lj = and e = — c 13 C. Acting on any velocity u M and any 
scalar ip (not necessarily the inflaton), the same gauge 
transformation leads to the matter perturbations 

Sip — ci^ip'C, Su = —ci^aC. (D10) 

Since by assumption the aether does not couple to mat- 
ter, and for the same reasons as in the adiabatic case, 
we expect Eqs. (D9) and (D10) then to be a solution 
of the matter equations of motion, no matter what the 
aether perturbation C actually is. Of course, for arbi- 
trary values of C, we cannot expect the ansatz to satisfy 
Einstein's equations, since the aether does couple to grav- 
ity. Inspection of the latter however reveals that the %, 
°i and diagonal 1 ; equations only contain spatial gradi- 
ents of the aether field, which can be neglected in the 
long- wavelength limit. The only equation in which the 
aether perturbation is not negligible at long wavelengths 
is ( B4 ) , which is actually satisfied by the ansatz ( D9 1 . 



Hence, it only remains to find out what the aether per- 
turbation C is. Substituting Eq. (D9| into the aether 



field equation (B2) results in a differential equation for 



the yet undetermined aether perturbation, 



C" + 2HC 



a 



H z 



1 - 



C14 



H' 



C = 0. 
(DH) 
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This equation has two independent solutions, which when 



plugged into (D9) and (D10) give the to two independent 



isocurvature modes, for which ( = to = 0. None of these 
solutions can be adiabatic, as the adiabatic mode has 
e = C, while along these solutions e = — C13C (recall 
that C13 = — 1 is a singular case.) A measure of the 
non-adiabaticity of these modes is the difference in the 
e-folding number between surfaces comoving with aether, 
and those comoving with matter, which equals 

SN= (l + c 13 )HC = -i±^V, (D12) 

Cl3 

and thus differs from zero if C13 ^ — 1. Since along these 
solutions all matter components (aside from the aether) 
share the same velocity, the two modes describe a matter- 
aether isocurvature perturbation, which is the only kind 
of isocurvature perturbation that can be generated if the 
aether does not couple to matter. For a constant equation 
of state, these two isocurvature modes reproduce those 
found in Subsection IB 2 bl 

Note that this method of generating solutions would 
break down if the anisotropic stress of matter on large 
scales were not negligible, as would happen for instance 
if the matter sector contained a second aether field. 



Appendix E: CMB anisotropies in the vector sector 

The effect of vector perturbations on the amplitude 
of CMB anisotropies is easily estimated in the approxi- 
mation of a sharp transition between thermal equilibrium 
and complete transparency at the moment of decoupling. 
Before the transition, photons and baryons are approxi- 
mated as a perfect fluid, whereas after the transition the 
radiation will be described in terms of a distribution of 
free photons. 

The number of photons in a phase space cell can be 
written as 



dn = n(x, p) J^rfx fe TTdpi, 



(El) 



where x k are space coordinates and pk are the spatial 
components of the momentum. For a gas of free photons, 
the number density in phase space obeys the collisionless 
Boltzmann equation 



dn ^ dn dx k dn dpk 
dn dx k drj dpk dr\ 



0. 



(E2) 



Further, we assume that the distribution of photons trav- 
eling in a given direction 1 at any given point has the 
Planckian spectrum, 



n = n(E/T). 



Here, 



E 



-p^ = -a 1 p 



(E3) 



(E4) 



is the energy of a photon as measured by an observer 
at rest in the coordinates x. The four-velocity of this 
observer is given by w M = (— ffoo) _1 ^ 2 ^o ' anc ^ m ^ ne ^ as * 
equality we have used that —goo = a 2 is unperturbed in 
the linearized vector sector. The local temperature 



T = T (r 1 )+5T(r ] ,x,l) 



(E5) 



depends not only on position, but also on the direction 
of arrival of the photons, 



h=Pi/p, where p = {5 tj PiPj) 1/2 . 



(E6) 



Before decoupling, when the system is in thermal equi- 
librium, the temperature anisotropy is just a dipole, cor- 
responding to the local motion of the photon fluid. This 
is characterized by the four- velocity <5u M . Note that n is 
a scalar, and so T is defined in such a way that the ratio 
y = E/T transforms as a scalar. In the co-moving frame, 
where the fluid is at rest, we have 



Ec 



+ Su^pn E-6u i p l 



t,. 



(E7) 



where the co-moving temperature T c = To + §o(t), x) is 
isotropic, and we have used Su° — (to linear ord er in 
Su 1 ). Since y = E/T = E c /T c , it follows from (JET]) that 
at the time of decoupling 



ST 



(?7doc,x, 1) = 5 + a Su l k 



(E8) 



Later, after decoupling, the photons arriving from differ- 
ent directions at a given spacetime point have originated 
at different locations on the surface of last scattering, 
which leads to anisotropies also in the higher multipoles. 

The monopole and dipole components in ( E8 ) are re- 
lated to the perturbations in Tq and Tf, which can be 
obtained from the expression 



-fJ 



n(y) 



FX 
p o 



d 3 p. 



(E9) 



Here p stands for the spatial components of the momen- 
tum, with lower indices. Let us consider the perturba- 
tion in the energy density. This will be related to the 
monopole component in the temperature anisotropy. For 
vector perturbations, 



2 jr 2 

dgoi = a 



= S l 



2 6 9ij = a^g- n-- ' ■ /•";. 



(E10) 
(Ell) 



the linearized metric determinant is y/—g — a 4 , and the 
condition p M p M = leads to 



Po = -p(l - SHi - F 1 ' 3 hi,), 
The energy density of photons is given by 
1 



Per. 



_ rpO _ 



n(y)p p dpd 1. 



(E12) 



(E13) 
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We can now use that po — —aT y(l + 5T/T) and p = 
aT y(l + S% + F^ kl + ST/T) to eliminate p and p in 
favour of y. After simple manipulations, one obtains 



Pc 13 



P (0) 



1 



d 2 \8T 



(E14) 



Using (E12| and (E18) in (E17l we have 



(E19) 



Since vector perturbations do not change the energy den- 
sity, we have Sp Cl3 = 0. Therefore, using (E8) in (E14) we 
find Sq = 0. Hence, the temperature anisotropy (E8) for 



vector perturbations in the perfect fluid is purely dipolar: 

^(7 7dcc ,x,l) = -5 i Z i + ^i k. (E15) 

For later convenience, here we have expressed the result 
in terms of the velocity perturbation with lower indices, 
which is gauge-invariant. 

The evolution of the temperature anisotropy after de- 
coupling can be inferred from the Boltzmann equation. 
Defining 



E = p/a, 



(E16) 



we have d v (E /T Q ) = d x * {E Q /T ) = 0, and d Pk {E /T ) = 
(l k /p)(E /T ). Substituting (|E3]) in (|E2]), and linearizing 
in perturbations, it is straightforward to show that 



Here, d/di] = d v + kd x i is the total derivative along 
the line of sight, and primes indicate partial derivatives 
with respect to r\. The result is expressed in terms of 
the gauge-invariant c ombi nations (5T/T ) + F' • 1 and 
qi = p% _ g% (E19) can be integrated along the 

trajectory of the photons x(?]) = (r/ — 770)1, from the time 
of decoupling ?7dec to the present time 770, to obtain the 
temperature anisotropy which is observed at present: 



ST 



Sui 



dec 



F' ■ 1 ) = ( Q% + — k ) + I drj Q' • 1. 

(E20) 

Here we have used the initial condition determined by 
(E15). As we mention in Subsection VIA a non- 



vanishing velocity perturbation Sui cannot be generated 
as long as the perfect fluid description is valid, so we shall 
ignore Sui in Eq. ( E20 1 , and simply write 



drj k dx k 



SE 
Eo~ 



ST 



h dp k 
p dr] 



= 0, 



(E17) 



where SE = E—E = -~(po+p)/a. The geodesic equation 
reads 



dpk = 1 dg i 
drj 2p° dx k 



k Pi 



F l 



PiPj 



(E18) 



ST 
To~ 



r'lo 

= (Q • l)dec + / dvQ'-l. 



(E21) 



Here we have also dropped the dipole term at the time 
of observation, since this is always subtracted. 
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